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Abstract 


This  report  summarizes  the  work  completed  in  the  months  of  June  1996-December  1999  for 
the  principal  investigator’s  Young  Investigator  Program  grant.  As  such  the  report  represents 
a  final  summary  of  work  and  is  organized  as  follows.  First,  in  a  brief  summary  chapter,  a 
review  of  the  accomplishments  of  the  work  supported  from  June  1996  to  December  1999  is 
given  providing  an  overview  of  the  successes  during  that  period.  Next,  two  chapters  giving 
details  of  new  work  during  the  six  months  since  the  last  annual  report  are  presented.  In  these 
results  an  experimentally  useful  solution  for  the  stress  intensity  factor  history  at  the  corner 
of  an  impacting  punch  is  derived  for  isotropic  and  orthotropic  materials.  Then  that  solution 
is  compared  to  an  experiment  using  interferometry  to  measure  the  stress  intensity  factor 
history  at  the  corner  of  an  impacting  punch  in  both  isotropic  and  anisotropic  materials.  Good 
agreement  between  the  model  and  the  experiments  is  seen.  Consequently,  it  is  concluded 
that  the  solutions  derived  throughout  the  duration  of  this  work  are  accurate — mathematically 
speaking  the  fundamental  solutions  are  actually  exact — and  describe  the  loading  under  blunt 
impact  quite  well.  With  such  accurate  descriptions  of  the  stress  in  the  material  under  impact, 
future  work  can  now  focus  on  understanding  the  failure  in  the  material  under  these  same 


conditions. 
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Chapter  1 

Summary  of  Accomplishments 


An  analytical  approach  to  blunt  impact  problems  on  isotropic  and  anisotropic  materials  us¬ 
ing  dynamic  fracture  mechanics  was  shown  to  be  valid  over  the  duration  of  the  grant.  This 
approach  used  the  mathematical  and  experimental  techniques  usually  associated  with  dy¬ 
namic  fracture  mechanics  to  examine  impact  of  a  blunt  object  on  a  wide  class  of  materials. 
For  simplicity  sake,  situations  were  limited  to  two-dimensional  loading  and  geometries.  Al¬ 
though  this  offers  some  limitations  oh  the  results,  it  also  allows  both  the  derivation  of  exact 
solutions  describing  events  occurring  on  the  interior  of  impacted  materials  and  the  experi¬ 
mental  observation  of  those  events,  neither  of  which  is  readily  achieved  for  three-dimensional 
problems.  In  addition,  the  two-dimensional  results  offer  useful  insight  into  three-dimensional 
impact  events  and  provide  valuable  intuition  about  such  events. 

1.1  Previously  Reported  Work 

First,  through  studies  of  impact  on  metals,  it  was  shown  in  this  work  that  fracture  toughness 
can  be  used  to  characterize  failure  of  metals  by  shear  localization  under  impact  conditions. 
Simple  solutions  for  the  stress  intensity  factor  in  the  dynamic  punch  test  were  derived  using 
known  fundamental  solutions  from  dynamic  fracture  mechanics.  The  results  of  these  models 
were  published  as 

K.M.  Roessig  and  J.J.  Mason,  “Dynamic  Stress  Intensity  Factors  in  a  Two  Di¬ 
mensional  Punch  Test,”  Engineering  Fracture  Mechanics,  60,  No.  4,  pp  421-435, 
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1998 


The  conclusion  from  the  models  was  that  the  stress  intensity  history  around  the  corner  of 
an  impacting  punch  was  more  severe  than  that  around  a  notch  tip  in  a  Kalthoff  test.  This 
result  was  significant  in  light  of  the  fact  that  much  attention  has  been  directed  toward  the 
Kalthoff  test  as  a  measure  of  resistance  to  shear  localization  during  impact.  As  it  was  shown 
by  the  model,  actual  impact  can  be  worse  than  the  Kalthoff  test. 

The  solution  for  the  stress  intensity  factor  history  around  the  corner  of  a  punch  impacting 
on  an  isotropic  material  was  then  used  to  test  many  types  of  metals.  It  was  seen  that  shear 
band  initiation  at  the  corner  of  a  punch  was  quite  similar  to  that  seen  at  a  notch  tip  and 
that  initiation  occurred  at  the  same  value  of  stress  intensity  factor  in  both  cases.  These 
results  were  quite  surprising  but  none-the-less  lead  one  to  believe  that  a  fracture  mechanics 
approach  to  modeling  shear  localization  may  be  possible.  The  results  were  published  as 

K.M.  Roessig  and  J.J.  Mason,  “Adiabatic  Shear  Localization  in  the  Impact  of 
Edge  Notched  Specimens,”  Experimental  Mechanics,  38,  No.  3,  pp.  196-203, 

1998 

After  examining  isotropic  materials,  attention  was  turned  to  anisotropic  materials  and 
solutions  for  stress  intensity  factors  in  two  dimensional  orthotropic  materials  were  derived. 
Initially,  solutions  using  numerical  methods  were  found.  This  approach  involved  reducing  the 
governing  equations  (in  this  case  the  conservation  of  momentum  in  terms  of  displacement) 
to  a  Fredholm  integral  equation  for  the  Laplace  transform  of  the  stress  intensity  factor  his¬ 
tory.  Although  numerical  inversion  of  the  Laplace  transform,  found  by  solving  the  Fredholm 
equation  numerically,  can  be  tricky,  it  was  successfully  performed  for  the  case  of  normal 
and  shear  point  loads  applied  to  a  finite  crack.  These  results  were  published  separately;  the 
results  for  normal  loads  appeared  as 

C.  Rubio-Gonzalez  and  J.J.  Mason,  “Green’s  Functions  for  the  Stress  Intensity 
Factor  Evolution  in  Finite  Cracks  in  Orthotropic  Materials,”  to  appear  Int.  J. 
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Fracture,  1999 


and  the  results  for  shear  point  loads  appeared  as 

C.  Rubio-Gonzalez  and  J.J.  Mason,  “Response  of  Finite  Cracks  in  Orthotropic 
Materials  due  to  Concentrated  Impact  Shear  Loads,”  J.  Applied  Mechanics,  66, 
no.  2,  pp.  485-491,1999 

Once  the  problems  for  finite  cracks  had  been  solved  it  was  discovered  that  instead  of 
reducing  the  governing  equations  to  a  Fredholm  integral  equation,  one  could  reduce  them 
to  a  Weiner-Hopf  equation.  This  is  quite  an  advance  in  the  analysis  techniques  for  dynamic 
fracture.  When  studying  isotropic  materials,  it  is  common  to  reduce  important,  fundamental 
problems  to  a  Weiner-Hopf  equation  by  using  Helmholtz  displacement  potentials.  Using  a 
similar  technique  on  anisotropic  materials  was  not  possible  because  proper  Helmholtz  po¬ 
tentials  do  not  exist  for  anisotropic  materials;  they  only  exist  for  isotropic  materials.  Here, 
however,  we  were  able  to  reduce  the  governing  equations  for  displacements  to  a  Weiner-Hopf 
equation  directly,  without  the  use  of  displacement  potentials.  Hence,  a  whole  new  class  of 
fundamental  problems  in  dynamic  fracture  mechanics  could  be  solved.  We  published  a  sum¬ 
mary  of  our  method,  as  it  is  applied  to  uniformly  loaded,  semi-infinite  cracks  in  orthotropic 
materials,  in  the  prestigious  Journal  of  Mechanics  and  Physics  of  Solidsas 

C.  Rubio-Gonzalez  and  J.J.  Mason,  “Closed  Form  Solutions  for  the  Dynamic 
Stress  Intensity  Factor  at  the  Tip  of  Uniformly  Loaded  Semi-infinite  Cracks  in 
Orthotropic  Materials,”  to  appear  J.  Mechanics  and  Physics  of  Solids,  1998 

Then  we  looked  at  point  loads  applied  to  semi-infinite  cracks  in  orthotropic  materials.  These 
closed  form  solutions  serve  as  Green’s  functions  for  solutions  to  a  wide  class  of  problems  and 
differ  from  our  earlier  numerical  solutions  in  that  they  are  closed  form  and  apply  to  semi¬ 
infinite  cracks  rather  than  finite  cracks.  The  new  method  of  reducing  the  problem  to  a 
Weiner-Hopf  equations  was  used  to  find  solutions  for  normal  point  loads,  published  as  as 
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C.  Rubio-Gonzalez  and  J.J.  Mason,  “Dynamic  Stress  Intensity  Factor  Due  to 
Concentrated  Normal  Loads  on  Semi-infinite  Cracks  in  Orthotropic  Materials,” 
to  appear  Journal  of  Composite  Materials,  1999 

and  shear  point  loads,  published  as 

C.  Wang,  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Dynamic  Stress  Intensity  Factor 
on  Semi-infinite  Cracks  in  Orthotropic  Materials  Due  to  Concentrated  Shear 
Impact  Loads,”  submitted  to  Int.  J.  of  Solids  and  Structures,  1998. 

These  solution,  in  contrast  to  the  solutions  described  in  the  preceding  paragraph,  are  closed 
form  and  not  nearly  as  tricky  to  evaluate  numerically.  Next,  we  turned  our  attention  to 
propagating  cracks.  The  same  method  allowed  us  to  find  two  approximations  which  describe 
the  variation  of  the  stress  intensity  factor  as  a  function  of  crack  velocity,  one  for  isotropic 
materials  that  agreed  with  known  solutions  and  one  for  highly  orthotropic  materials  that 
contributed  new  knowledge  in  the  area  of  dynamic  fracture.  This  work  was  published  as 

C.  Rubio-Gonzalez  and  J.J.  Mason,  “Dynamic  Stress  Intensity  Factor  for  a  Prop¬ 
agating  Semi-Infinite  Crack  in  Orthotropic  Materials,”  to  appear  Int.  J.  of  En¬ 
gineering  Science,  1998. 

Such  results  are  useful  in  characterizing  material  resistance  to  dynamic  crack  propagation. 

Having  found  a  new  mathematical  technique  for  solving  two  dimensional  dynamic  frac¬ 
ture  problems  in  orthotropic  materials,  we  attempted  to  apply  the  technique  to  quasi-three- 
dimensional  problems.  We  examined  the  loading  of  penny  shaped  cracks  in  transversely 
isotropic  materials.  Using  an  approach  based  on  the  approach  derived  for  two-dimensional 
problems,  we  were  able  to  approximate  the  penny  shaped  crack  problem  by  a  simpler  problem 
that  could  be  reduced  to  a  system  of  two  Weiner-Hopf  equations.  This  problem  was  solved 
and  the  results  show  how  a  penny  shaped  crack  might  grow  if  it  were  loaded  dynamically  in 
tension  or  shear.  The  paper  resulting  from  the  analysis  is  published  as 
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C.Y.  Wang  and  J.J.  Mason,  “The  Dynamic  Stress  Intensity  Factor  and  Strain 
Energy  Release  Rate  for  a  Semi-infinite  Crack  in  Rotated  Transversely  Isotropic 
Materials  due  to  Uniform  Impact  Loading,”  submitted  to  Int.  J.  of  Fracture, 

1999. 

This  results  is  useful  when  studying  the  initiation  of  delaminations  in  composites  under 
impact  loading. 

1.2  Recent  Work 

Finally,  during  the  last  six  months  of  the  funded  effort,  experimental  investigation  of  the 
application  of  the  solutions  for  stress  intensity  factors  in  orthotropic  materials  to  dynamic 
crack  initiation  in  composites  under  impact  conditions  was  completed.  First,  the  fundamen¬ 
tal  solutions  derived  for  semi-infinite  cracks  in  orthotropic  materials  were  used  to  calculate 
the  stress  intensity  factor  in  the  dynamic  punch  test.  This  work  is  presented  in  Chapter  2 
and  has  been  submitted  as 

C.  Rubio-Gonzalez  and  J.J.  Mason,  “Elastodynamic  Analysis  of  Finite  Punch 
and  Finite  Crack  Problems  in  Orthotropic  Materials,”  submitted  to  Theoretical 
and  Applied  Fracture  Mechanics,  1999 

With  that  solution  in  hand,  we  performed  experiments  on  isotropic  and  orthotropic  materials 
and  measured  the  stress  intensity  factor  history  in  the  dynamic  punch  tests.  This  work  is 
presented  in  Chapter  3  and  published  as 

C.  Rubio-Gonzalez  and  J.J.  Mason,  “Experimental  Investigation  of  Dynamic 
Punch  Tests  on  Isotropic  and  Composite  Materials,”  submitted  to  Experimental 
Mechanics,  1999 

The  results  agree  quite  well  with  the  solution  in  Chapter  2  and  thus  verify  solutions  presented 
here  and  the  method  used  to  derive  them. 
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1.3  Summary  of  Journal  Publications 

Overall,  the  work  has  been  quite  successful  and  has  lead  to  significant  advances  in  the 
modeling  of  dynamic  fracture  in  anisotropic  materials.  A  total  of  eleven  journal  publications 
have  been  produced  and  fifteen  conference  presentations  and  invited  lectures  have  been 
given  by  the  principal  investigator.  Also,  two  students  have  received  their  Ph.D.  degrees; 
one  from  an  under-represented  group  in  engineering.  Future  work  should  be  directed  toward 
using  these  solutions  in  experimental  investigation  of  material  resistance  to  fracture  under 
dynamic  loading  conditions  such  as  underwater  blasts  or  fragment  impact. 

The  following  list  gives  the  journal  publications  resulting  from  the  funded  work. 

1.  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Elastodynamic  Analysis  of  Finite  Punch  and 
Finite  Crack  Problems  in  Orthotropic  Materials,”  submitted  to  Theoretical  and  Applied 
Fracture  Mechanics,  1999 

2.  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Experimental  Investigation  of  Dynamic  Punch 
Tests  on  Isotropic  and  Composite  Materials,”  submitted  to  Experimental  Mechanics, 
1999 

3.  C.Y.  Wang  and  J.J.  Mason,  “The  Dynamic  Stress  Intensity  Factor  and  Strain  Energy 
Release  Rate  for  a  Semi-infinite  Crack  in  Rotated  Transversely  Isotropic  Materials  due 
to  Uniform  Impact  Loading,”  submitted  to  Int.  J.  of  Fracture,  1999 

4.  C.  Wang,  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Dynamic  Stress  Intensity  Factor  on 
Semi-infinite  Cracks  in  Orthotropic  Materials  Due  to  Concentrated  Shear  Impact  Loads,” 
submitted  to  Int.  J.  of  Solids  and  Structures,  1 998 

5.  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Dynamic  Stress  Intensity  Factor  for  a  Propagating 
Semi-Infinite  Crack  in  Orthotropic  Materials,”  to  appear  Int.  J.  of  Engineering  Science, 
1998 
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6.  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Dynamic  Stress  Intensity  Factor  Due  to  Con¬ 
centrated  Normal  Loads  on  Semi-infinite  Cracks  in  Orthotropic  Materials,”  to  appear 
Journal  of  Composite  Materials,  1998 

7.  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Closed  Form  Solutions  for  the  Dynamic  Stress 
Intensity  Factor  at  the  Tip  of  Uniformly  Loaded  Semi-infinite  Cracks  in  Orthotropic 
Materials,”  to  appear  J.  Mechanics  and  Physics  of  Solids,  1998 

8.  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Green’s  Functions  for  the  Stress  Intensity  Factor 
Evolution  in  Finite  Cracks  in  Orthotropic  Materials,”  to  appear  Int.  J.  Fracture,  1999 

9.  C.  Rubio-Gonzalez  and  J.J.  Mason,  “Response  of  Finite  Cracks  in  Orthotropic  Mate¬ 
rials  due  to  Concentrated  Impact  Shear  Loads,”  J.  Applied  Mechanics,  66,  no.  2,  pp. 
485-491,1999 

10.  K.M.  Roessig  and  J.J.  Mason,  “Adiabatic  Shear  Localization  in  the  Impact  of  Edge 
Notched  Specimens,”  Experimental  Mechanics,  38,  No.  3,  pp.  196-203,  1998 

11.  K.M.  Roessig  and  J.J.  Mason,  “Dynamic  Stress  Intensity  Factors  in  a  Two  Dimensional 
Punch  Test,”  Engineering  Fracture  Mechanics,  60,  No.  4,  pp  421-435,  1998 
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Chapter  2 

Elastodynamic  Analysis  of  Finite 
Punch  and  Finite  Crack  Problems  in 
Orthotropic  Materials 
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Overview 


The  transient  elastodynamic  response  of  the  finite  punch  and  finite  crack  problems  in  or¬ 
thotropic  materials  is  examined.  Solution  for  the  stress  intensity  factor  history  around  the 
punch  corner  and  crack  tip  is  found.  Laplace  and  Fourier  transforms  together  with  the 
Wiener-Hopf  technique  are  employed  to  solve  the  equations  of  motion  in  terms  of  displace¬ 
ments.  A  detailed  analysis  is  made  in  the  simplified  case  when  a  flat  rigid  punch  indents 
an  elastic  orthotropic  half-plane,  the  punch  approaches  with  a  constant  velocity  normally 
to  the  boundary  of  the  half-plane.  An  asymptotic  expression  for  the  singular  stress  near 
the  punch  corner  is  analyzed  leading  to  an  explicit  expression  for  the  dynamic  stress  inten¬ 
sity  factor  which  is  valid  for  the  time  the  dilatational  wave  takes  to  travel  twice  the  punch 
width.  In  the  crack  problem,  a  finite  crack  is  considered  in  an  infinite  orthotropic  plane.  The 
crack  faces  are  loaded  by  impact  uniform  pressure  in  mode  I.  An  expression  for  the  dynamic 
stress  intensity  factor  is  found  which  is  valid  while  the  dilatational  wave  travels  the  crack 
length  twice.  Results  for  orthotropic  materials  are  shown  to  converge  to  known  solutions  for 
isotropic  materials  derived  independently. 

2.1  Introduction 

The  punch  problem  is  of  great  importance  in  solid  mechanics  for  its  multiple  technical  appli¬ 
cations  including  ballistic  impact,  explosives,  metal  forming  and  manufacturing  operations 
such  as  punching  and  blanking.  While  the  quasi-static  punch  problem  is  a  well  studied  field 
in  contact  mechanics  [1,2],  there  is  little  work  on  the  dynamic  case.  In  the  dynamic  case  the 
punch  approaches  the  material  with  a  certain  velocity,  hence  wave  propagation  is  involved 
complicating  the  mathematical  analysis. 

Brock  [3,  4]  considered  the  problem  of  a  rapid  indentation  of  an  isotropic  half-plane  by  a 
smooth,  flat  and  rigid  semi-infinite  punch.  The  high  stresses  generated  under  the  punch  edge 
were  assumed  to  generate  slip  mechanisms,  and  the  zone  was  allowed  to  relax  the  stresses 
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under  the  punch  edge.  That  is,  the  stress  relaxation  feature  of  the  zone  was  imposed  to 
remove  the  square-root  singularity  in  the  normal  contact  stress  under  the  punch  edge,  and 
this  condition  led  to  simplify  the  mathematical  treatment.  It  has  been  shown,  however,  in 
a  series  of  experiments  conducted  in  isotropic  and  unidirectional  composite  materials  [5], 
that  such  plastic  zone  is  small  enough  like  that  usually  assumed  near  a  crack  tip  in  fracture 
mechanics,  and  the  singular  stress  field  prevails  near  the  punch  edge.  In  those  experiments,  a 
rigid  flat  projectile  with  sharp  corners  impacted  (at  moderate  velocities  ~20  m/s)  one  side  of 
rectangular  specimens,  and  the  stress  intensity  factor  time  evolution,  which  is  a  parameter 
related  to  the  singular  normal  stress  near  the  punch  edge,  was  measured  using  the  laser 
interferometry  method  of  coherent  gradient  sensor  (CGS  method),  in  combination  with  high 
speed  photography. 

Punch  geometries  other  than  flat,  have  been  considered  even  for  general  anisotropic 
materials,  e.g.  parabolic  and  wedge  shaped  punches,  [6,  7].  In  those  problems  the  contact 
area  grows  at  constant  speed  and  the  self-similar  feature  of  the  dynamic  fields  is  exploited 
using  the  general  methodology  derived  by  Willis  [6]  for  such  kind  of  problems. 

The  dynamic,  flat  finite  punch  problem  of  width  21  in  isotropic  materials  was  considered 
by  Roessig  and  Mason  [8].  They  solved  an  equivalent  problem  where  the  rigid  punch  is  re¬ 
placed  by  a  compressive  wave  impinging  two  semi-infinite  external  cracks.  The  computation 
of  the  stress  intensity  factor  history  for  the  time  interval  0  <  t  <  Al/cd,  was  carried  out  by 
integration  of  a  Green’s  function  in  space  and  time.  The  appropriate  Green’s  function  is  the 
solution  derived  by  Freund  [9]  for  the  semi-infinite  crack  problem  with  impact  point  loads  ap¬ 
plied  normally  on  the  crack  faces  at  certain  distance  Ip  from  the  crack  tip.  Even  though  this 
procedure  may  be  extended  for  the  case  of  orthotropic  materials,  since  the  needed  Green’s 
function  has  been  developed  by  Rubio-Gonzalez  and  Mason  [10],  another  approach  will  be 
used  in  this  work.  A  procedure  similar  to  that  presented  by  Poruchikov  [11]  is  employed 
here  which  allows  the  determination  of  an  explicit  expression  for  the  stress  intensity  factor 
history  in  the  time  interval  0  <  t  <  Al/cd- 
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In  the  following,  the  semi-infinite  and  finite  punch  problems  are  considered.  A  more 
general  semi-infinite  punch  problem  is  solved  first  where  a  non-homogeneous  boundary  con¬ 
dition  outside  the  punch  is  allowed.  This  strategy  will  permit  to  construct  a  solution  for  the 
finite  punch  problem  as  a  superposition  of  semi-infinite  punch  solutions.  The  expression  for 
the  dynamic  stress  intensity  factor  in  the  finite  punch  problem  in  orthotropic  materials  is 
new  in  the  literature.  The  dynamic  stress  intensity  factor  for  the  semi-infinite  punch  coin¬ 
cides  with  that  for  the  semi-infinite  crack  under  uniform  impact  loads  on  the  crack  faces  in 
mode  I  and  it  was  derived  first  in  [12]. 

The  same  strategy  is  applied  then  to  construct  a  solution  for  the  finite  crack  problem 
in  orthotropic  materials.  The  stress  intensity  factor  history  Ki{t)  is  determined  for  the 
time  interval  0  <  t  <  Al/cd,  being  21  the  crack  length.  The  dynamic  crack  problem  in 
orthotropic  materials  was  considered  first  by  Bandyopahya  and  Kassir  [13]  using  a  method 
derived  by  Chen  and  Sih  [14].  Such  method  gives  an  approximate  solution  for  Ki{t)  after 
performing  a  numerical  solution  of  a  Fredholm  integral  equation  and  numerical  Laplace 
transform  inversion.  Nevertheless  this  method  permits  the  determination  of  Ki{t)  for  long 
times,  it  may  be  inaccurate  if  the  proper  numerical  schemes  are  not  selected,  as  will  be 
illustrated  in  this  work  comparing  Kassir  and  Banyopadhyay  [13]  results  with  the  exact 
solution  developed  here.  The  Wiener-Hopf  technique  has  been  applied  in  the  finite  crack 
problems  in  isotropic  materials  [15, 16].  Freund  [16]  developed  an  expression  for  the  dynamic 
stress  intensity  factor  using  a  fundamental  solution  for  a  moving  dislocation  which  is  a 
different  approach  to  that  used  in  this  work.  It  is  noted  that  Freund’s  solution  for  the 
isotropic  case  is  recovered  from  the  orthotropic  expression  developed  here  with  the  proper 
substitution  of  the  elastic  constants. 
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2.2  Semi-infinite  Punch 

2.2.1  Governing  Equations 

Consider  that  at  t  =  0  a  smooth  punch  begins  its  indentation  into  an  elastic  orthotropic 
half-plane  y  >  0,  at  a  semi-infinite  interval  x  >  0  of  its  boundary,  see  figure  2.1.  Let  Eu  t^ij 
and  Uij  {i,j  =  1,2,3)  be  the  engineering  elastic  constants  of  the  material  where  the  indices 
1,  2,  and  3  correspond  to  the  directions  (x,  y,  z)  of  a  system  of  Cartesian  coordinates  chosen 
to  coincide  with  the  axes  of  material  orthotropy. 


Figure  2.1:  Schematic  of  the  semi-infinite  punch. 


The  problem  is  restricted  to  two  dimensions  with  wave  propagation  in  the  x  -  y  plane 
only.  By  setting  all  the  derivatives  with  respect  to  2:  to  be  zero,  it  is  readily  shown  that  the 


displacement  equations  of  motion  [17]  reduce  to 

d'^u  _  1  d‘^u 

d^v  ,,  ,  d^u  1  d‘^v 


(2.1) 

(2.2) 


where  u  and  v  are  the  x  and  y  components  of  the  displacement  vector  and  Cn,  C12  and  C22 


are  non-dimensional  parameters  related  to  the  elastic  constants  by  the  relations. 


El 

“  M12[1  -  (E2/El)*'f2l’ 

C22  =  (.E2/E,)cu,  (2-3) 

C12  =  1^12^22  =  ^'2lCll, 
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for  generalized  plane  stress,  and  by 


Cll 

C22 

C12 

A 


(1  - 1^231^32), 
(1  — 1^131^31)) 


El 

fJ>12^ 

E2 

fll2^ 

El  ,  ^E2  . 

1  —  l'i2l'21  ~  ^23^32  ~  ^^31*^13  ~  ^'12^^23^31  ~~  ^13^21^32 1 


(2.4) 


for  plane  strain.  In  the  orthotropic  solid,  Cg  =  sjliyil P  represents  the  velocity  of  the  in-plane 
shear  wave  propagating  along  the  the  principal  material  axes  and  p  is  the  mass  density.  The 
stresses  are  related  to  the  displacements  by  the  equations: 


O’! 

du  dv 

Pl2 

=  Cll- - f-Ci2^, 

ox  ay 

(Jy 

du  dv 

y 

P‘12 

=  Ci2-^ - I-C22-5-) 

dx  dy 

'^xy 

du  dv 
dy  dx 

Pl2 

(2.5) 


2.2.2  Method  of  Solution 

Let  /o(a:,  t)  be  a  function  describing  the  motion  of  the  punch.  Thus  the  boundary  conditions 
are 


ay{x,0,t)  =  0  for  -  oo  <  a:  <  0, 

Txy{x,0,t)  =  0  for  -  oo  <  a;  <  oo,  (2.6) 

u(a;,0,t)  =  fo{x,t)  for  a;  >  0. 

In  addition,  the  condition  of  zero  displacements  at  infinity  and  zero  initial  conditions  are 
assumed.  Note  that  the  first  two  equations  in  (2.6)  indicate  that  the  boundary  y  =  0  outside 
the  punch  is  stress  free.  To  ensure  the  uniqueness  of  the  solution  [11]  it  is  necessary  to 
restrict  the  behavior  of  the  displacement  in  the  vicinity  of  the  edge  of  the  punch  (a:  =  0), 
that  is, 

v{x,y,t)  =  v{t)  +  0{r^),  e  >  0,  r 0,  {r  =  y/x^  +  y^).  (2.7) 
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Instead  of  solving  the  problem  with  the  homogeneous  boundary  condition  <jy{x,  0,  t)  =  0 
in  equations  (2.6),  consider  a  more  general  problem  with  a  given  non-homogeneous  boundary 
condition  ay{x,0,t)  =  a-{x,t)  for  rr  <  0,  i.e.,  equations  (2.6)  become 

ay{x,0,t)  =  a-{x,t)  for  -  oo  <  x  <  0, 

Txy{x,0,t)  =  0  for  -oo<x<oo,  (2.8) 

u(x,0,t)  =  for  a;  >  0, 

where  <t_(x,  t)  is  assumed  to  be  equal  to  zero  at  any  x  <  — —  k{k>  0).  The  new  problem 
with  boundary  conditions  (2.8)  will  be  useful  in  solving  dynamic  problems  with  the  finite 
punch.  Note  that  with  a-(x,t)  =  0  the  standard  problem  with  the  semi-infinite  punch  is 
recovered.  The  function  cr_(x,t)  is  limited  everywhere  except  for  the  point  x  =  0  where  it 
should  satisfy  the  constraint 

a-(x,t)  =  0(x"'^'^*),,  e  >  0,  x  — 0  .  (2-9) 


The  method  of  solution  of  the  governing  equations  presented  here  follows  that  described 
by  Poruchikov  [11]  for  the  isotropic  case  with  some  significant  differences.  Displacement 
potentials  are  not  used.  In  equations  (2.1)  and  (2.2),  the  time  variable  may  be  removed  by 
application  of  the  Laplace  transform 

/•(p)  =  1°°  m  dt,  Sit)  =  ^  /-(p)  e"  dt,  (2.10) 


where  Br  denotes  the  Bromwich  path  of  integration  which  is  a  line  parallel  to  the  imaginary 
axis  in  the  p-plane.  Applying  relations  (2.10)  to  equations  (2.1)  and  (2.2)  and  using  zero 
initial  conditions  for  the  displacements  and  velocities,  the  transformed  field  equations  become 


d'^u*  &^u*  .  d'^v* 


d^v*  d'^v* 
-rrr  +  C22-7rT 


+  (1  +  C12) 


dxdy 


-Kv 


(2.11) 

(2.12) 


where  the  transformed  displacement  components,  u*  and  v*,  are  now  functions  of  the  vari¬ 


ables  X,  y,  and  p.  The  application  of  the  Laplace  transform  to  the  boundary  conditions  (2.8) 
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gives 


o-*(a;,0,p)  =  <y*_{x,p)  for  -  oo  <  a;  <  0, 
r*y(x,0,j9)  =  0  for  -  oo  <  a:  <  oo,  (2.13) 

v*{x,0,p)  =  fo{x,p)  for  x>0. 


To  obtain  a  solution  of  the  differential  equations  (2.11)  and  (2.12)  subject  to  conditions 
(2.13),  the  Fourier  transform  is  applied, 

/oo  1  fOO 

fix)  dx,  fix)  =  —  Fis)  ds.  (2.14) 

-OO  J—oo 

It  is  assumed  that  the  displacements  in  the  Laplace  transform  domain  have  the  form 


1  ro® 

u*ix,y,p)  =  7;-/  Ais,y,p)e~^^^  ds,  (2.15) 

ZTT  ^—00 
1  roo 

v*ix,y,p)  =  —  /  jB(s,y,p)e““'"  ds,  (2.16) 

ZTT  J—oo 


where  A  and  B  are  the  Fourier  transforms  of  the  Laplace  transform  of  the  displacements, 
u*  and  V*,  respectively,  and  are  yet  to  be  determined.  Substituting  these  transforms  into 
equations  (2.11)  and  (2.12),  the  functions  A  and  B  are  found  to  satisfy  the  simultaneous 
ordinary  differential  equations 

icns"^ +p^/<?^A-^  +  il-\rCu)is^  =  0,  (2.17) 

^2  p  ^  A 

(s2+pVc2)5-C22-^  +  (l  +  Ci2)is—  =  0.  (2.18) 

The  solution  of  these  equations  which  vanishes  for  y  — 00  is 


Ais,y,p)  =  Aiis,p)e  +  A2is,p)e 
Bis,y,p)  =  — ^>li(s,p)e“^‘^  -  — ^2(s,p)e"'^2^, 


(2.19) 


where  Ai  and  A2  are  arbitrary  functions  and  ajis,p)  stands  for  the  functions 


(1  +  Cl2)7i 


(2.20) 
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with  7i  and  7!  being  two  distinct  roots  of  the  quadratic  equation 
C227^  +  [(ci2  +  2ci2  —  CiiC22)s^  -  (1  +  C22)j^ / +  (cnS^  +  / Cs)(s^  +  P^/c,)  =  0.  (2.21) 


The  branches  of  7i(p,s)  and  j2{p,s)  are  chosen  in  the  appropriate  Riemann  sheets 
and  having  positive  real  parts.  It  has  been  shown  by  Willis  (1973)  and  Poruchikov  (1993) 
that  such  a  restriction  on  71  and  72  is  possible  for  the  general  anisotropic  case  as  long  as 
Re{p}  >  0.  This  fact  is  a  consequence  of  the  positive  definiteness  of  the  strain  energy.  The 
expressions  for  the  displacements  in  the  Laplace  transform  domain  become: 

u*  =  —  r  +  A2e-'^^^)e-^‘’^  ds,  (2.22) 

27r  J-00 

_ A  rOO  ^—isX 

V*  =  —  {aiAie-^^y  +  a2A2e-''^y) - ds,  (2.23) 

27r  J-00  s 

and  using  (2.5)  the  corresponding  expression  for  T*y  is  given  by 

r*  =  r  [(ai  +  71)^16-^^"  +  {OC2  +  72)A2e-^^^]e-^^^  ds.  (2.24) 

^  Ztt  J —00 

Applying  the  second  condition  of  (2.13)  to  equation  (2.24)  yields 

A2{s,p)  =  -^iAi(s,p), 

=  (2.25) 

0(2  +  72 

Therefore,  the  expressions  for  the  transformed  components  of  displacement  become 

u*{x,y,p)  =  ^  r  -  Pie-'^^nA,{s,p)e-^^^  ds,  (2.26) 

Ztt  J  ~oo 

V(x,y,p)  =  ^  r  {ate-’'”  -  ds,  (2.27) 

ZTT  J-oo  S 

and  the  associated  stress  components  are  given  by 

a*  =  f  [(ciis^-ai7ici2)e"^'^-(cns^-a272ci2)/5ie”'^"^]^^^^j^e“*®"^c(:^.28) 

ZTT  J-oo  ^ 

<  =  r  [{cuS^  -  -  {cns'‘ - 

^  ZTT  J-oo  ^ 

T-  =  r  (a,+'„)\e-’'’'  -e-">]Ai{s,p)e-‘“ds.  (2.30) 

^  ZTT  J—oo 
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The  following  functions  are  introduced, 


E{s,p)  =  -(ai  -  l3ia2)Ai{s,p), 
s 

F{s,p)  =  -■ r - a - 4[C12S^-Q!i7iC22- A(Cl2S^-tt272C22)], 

(o-i  -  Pia2)^ 


V/^-2^  +  V- 

C11C22  ~  C12 
C22  ’ 


(2.31) 

(2.32) 

(2.33) 

(2.34) 


where  the  velocity  ca  =  ^/^Cs  represents  the  dilatational  wave  speed  along  the  x— axis.  In 
view  of  the  first  and  third  boundary  conditions  in  (2.13),  equation  (2.27)  and  (2.29)  now  yield 
the  following  pair  of  dual  integral  equations  for  the  determination  of  the  function  E{s,p) 

<fe0.p)  =  ^r  E{s,p)e->-‘ ds  =  al(x,p)  -  cx=  <  x  <  (5.35) 

27r  J-oo  ^s^-\-p^/(?^ 

_ 7  TOO 

v'{x,0,p)  =  ^  E{s,p)e-'‘‘ds  =  /;{x,p)  0  <  X  <  oo,  (2.36) 

Ztt  J—oo 

Let  v*_{x,p)  be  the  unknown  Laplace  transform  of  the  vertical  displacement  on  the 
negative  x— axis,  and  <7^(x,p)  be  the  unknown  Laplace  transform  of  the  normal  stress  on 
the  positive  x— axis,  so  that 

*/  n  \  f  /o(^)P)  for  X  >  0  */  n  \  f  a:  >  0  (n 

*’  (*'0.P)  =  |„l(x;p)  for  x<0,  for  X<0. 

then  we  can  write  the  Laplace  transform  of  the  normal  stress  and  the  vertical  displacement 

on  the  whole  boundary  y  =  0  as 

£;(s,p)e-"-  ds  =  a*_{x,p)  +  al{x,p),  (2.38) 

zTT  J —oo  ^  /  c2  _1_  r>2  //^2 


F{s,p) 

£'(s,p)e  """  ds  =  crl(x,p) +  cT+(x,p), 

(2.38) 

yjs^+p^/cl 

-i 

27r  J —oo 

E{s,p)e~^^^  ds  =  u*  (x,p)  + /o(x,p). 

(2.39) 

and  by  Fourier  transform  inversion,  these  equations  give 

ipi^i-jEkM=E(s,p)  =  E-(s)  +  E+(s). 

-iE{s,p)  =  V-{s)  +  V+{s), 


(2.40) 

(2.41) 


where 

E_(s)=  r  a*_{x,p)e^^^dx  E+(s)  =  1°°  aX{x,p)e^^^dx,  (2.42) 

J— CO  *'0 

V-{s)=  r  vl{x,p)e^^^dx  V+is)  =  r  f^{x,p)e^^^dx,  (2.43) 

J—oo  *'0 

Prom  a  priori  knowledge  of  properties  of  solutions  of  elastic  wave  equations  it  is  possible 
to  show  that  cr;^(a:,p)  and  v^{x,p)  are  exponentially  bounded  at  infinity  and  this  ensures 
the  existence  of  their  Fourier  transform  (2.42)  and  (2.43).  The  analysis  is  similar  to  that 
presented  by  Rubio-Gonzalez  and  Mason  [12]  for  the  crack  problem.  In  addition  it  is  shown 
by  Noble  [18]  that 

if  \aKx,p)\  <  as  a;  ->  +oo  then  I!+(s)  is  analytic  in  Im(s)  =  A  >  A_, 

and  if  |u!.(x,p)|  <  M2e'^+®  as  x  -oo  then  V_(s)  is  analytic  in  Im(s)  =  A  <  A+. 

Since  all  disturbances  are,  at  least  limited  when  x  ->  +oo  and  are  absent  if  x  <  —cat  —  k, 
then  <t;(x,p)  ~  0(1)  as  x  -4-  +oo  and  v*_{x,p)  ~  0(eP^/‘=‘')  as  x  -oo.  Therefore,  it  is 
possible  to  show  that  A_  =  0  and  A+  =  p/cd-  Thus,  the  Fourier  transforms  (2.42)  and  (2.43) 
lead  to  conclude  that  E+(s)  and  VI  (s)  will  be  analytic  in  Im(s)  >  0  and  Im(s)  <  p/cd 
respectively. 


2.2.3  Wiener-Hopf  technique 


Eliminating  E{s,p)  from  (2.40)  and  (2.41)  we  obtain  a  Wiener-Hopf  equation 

E_(s)  +  S+(s)  =  -W2«-7==7?  !''-(*)  +  V'+(s)l  (2-44) 

^s^+p^/cl 

which  contains  only  the  two  unknown  functions  E+(s)  and  VT(s),  and  now  the  Wiener-Hopf 
technique  can  be  applied  as  follows.  Suppose  that  the  function  Lis')  is  defined  and  factored 


as 

L(s) 

then  equation  (2.44)  becomes 


L-js)  ^  Fjs.p) 

^+i^) 


(2.45) 


E_(s)L+(s)  -t-  E+(s)L+(s)  =  -L_(s)y_(s)  -  L_(s)14(s).  (2.46) 
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Assume  that  the  function  D{s)  is  defined  and  decomposed  as 


D(s)  =  E-(s)i+(s)  +  L_(s)nW 


(2.47) 


=  D^{s)  +  D-{s), 


then  equation  (2.46)  becomes 


D+{s)  +  E+(s)L+(s)  =  -L.{s)V.is)  -  D.{s)  =  Py(s). 


(2.48) 


The  first  member  of  this  equation  is  analytic  in  the  upper  half  of  the  s-plane  Im(s)  >  0  and 
the  second  member  is  analytic  in  the  lower  half  of  the  s-plane  Im(s)  <  p/ca-  Therefore,  the 
regions  of  analyticity  overlap.  Using  the  Liouville’s  theorem  to  determine  W(s),  solutions 
for  E+(s)  and  U-(s)  can  be  found. 

After  algebraic  manipulation  it  is  found  that  the  function  F{s,p)  reduces  to 


F{s,p) 


+  pVcKti  +  72)(1  +  C12) 

(ciiS^  +  P^/cl)  [c22(7i  +  72)  +  7172(1  +  C12)C22  + 

S^(C12  +  C12  -  C11C22)  -  C22pVc^]  } 


{c22Ci27i72  ~  7172(1  +  Ci2)Ci2S^  + 


(2.49) 


The  only  zeros  of  F{s,p)  are  of  the  form  s  =  dzip/cR  where  cr  is  the  Rayleigh  wave  speed. 
This  can  be  seen  by  substituting  s  =  ip/v  in  F{s,p),  letting  F{s,p)  =  0  and  dividing  by  the 
non-zero  factors,  then  F{ip/v,p)  =  0  reduces  to 


( C11C22  -  cfa 


which  is  the  Rayleigh  function  for  orthotropic  materials  [19].  The  roots  of  this  function  are 
V  =  ±cr.  For  a  proof  that  equation  (2.49)  has  only  these  two  zeros  see  [12] 


Consequently,  the  first  step  in  factoring  L{s)  is  to  define 


F(s)  = 


F{s,p) 


(2.50) 


It  can  be  shown  that  .F(s)  — 1  as  s  — >■  00,  (the  constant  ^  in  (2.32)  was  chosen  to  make  this 
possible).  The  function  F{s)  is  regular  and  F(s)  7^  0  in  the  s— plane  cut  as  shown  in  figure 
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Im{s) 

* 


Figure  2.2:  Branch  cuts  of  F{s)  in  the  s— plane. 


2.2,  the  only  singularities  are  the  branch  points  shared  with  71  and  72.  Where  the  branch 
points  of  7i  and  72  are 


for  7i; 
for  72;  s 


ip 

s  =  ±— , 

^5 


±- 


=  ±^ 


It  is  well-known  that  factorization  is  accomplished  most  directly  for  functions  that  ap¬ 
proach  unity  as  |s|  ^  00  and  that  have  neither  zeros  nor  poles  in  the  finite  plane;  F{s)  is 
an  example  of  such  a  function.  Therefore,  using  Cauchy’s  integral  formula  it  can  be  shown 
that  [16] 

=  exp(^/ 

[  ZTTZ  Jr±  z  —  s  J 


/r±  z  -  s 

where  F{s)  =  F+(s)F_(s)  and  r_  (r+)  is  the  contour  enclosing  the  branch  cut  between 


+ip/cd  and  +ip/cs,  {-ipjcd  and  -ip/cs).  Using  the  fact  that  F{s)  —  F{s),  one  can  write 

^  f  — 1  /-i/Ci  ,  /Im[F(ipu;)]\  dw 

F±(s)  =  exp  ^  —  /  tan-^  - u 

I  TT  Ji/cd  VR'e[F(zpu;)]/  u;  T  ^ 

Note  that  by  making  s  =  ipC,  in  this  equation,  F±{s  =  ipC}  =  F±(C)  becomes  a  function  only 
of  C-  That  is 

(  — 1  r\/cjt  (jy) 

'  (2.51) 


r_i  fi/c^ 

F±(0  =  exp  {  —  /  tan  [R{w)] 

(  TT  Jl/Cd 


W±C 
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where 


Re[/'(ipiy)] 

Returning  to  the  factorization  of  L{s)  we  have 

T(^\  -  hz!^-n 

~  r  A*12S  r—  „  ,  « 


therefore 


L+(s) 


+  P^/4 


=  M12^ 


F+(s)F_(s)(s  +  ip/cR){s  -  ip/cR) 
Js  +  ip/cdJs  -  ip/cd 


(2.52) 


(2.53) 


L_(s)  = 


■^+(5)  — 


yJS  -  ip/Cd 

fs  +  ip/cd  1 


(2.54) 


(2.55) 


-  {s  +  ip/cR)FM 

Using  the  sum  splitting  formula  [18]  for  the  function  D{s)  defined  in  (2.47)  it  is  found  that 


D^{s) 


1  roo+iXo  D{z) 


''  ''  27ri  J-oo+iXo  z  —  s 
D.{s)  =  D{s)-D+{s) 


(2.56) 

(2.57) 


where  Aq  is  such  that  A_  <  Aq  <  A+. 

Each  side  of  equation  (2.48)  is  analytic  in  one  of  the  overlapping  half  planes,  and  the  sides 
coincide  on  the  strip  of  overlap.  Consequently  each  side  of  (2.48)  is  the  analytic  continuation 
of  the  other  into  its  complementary  half  plane;  so  that  the  two  sides  together  represent  one 
and  the  same  entire  function  W^(s).  The  entire  function  will  be  determined  by  its  behavior  at 
|s|  — >  00  which  is  related  with  the  behavior  of  physical  quantities  near  a;  =  0.  First  note  that 
L+(s)  ~  and  L-{s)  ~  as  |s|  ->  00.  In  addition,  it  may  be  proven,  according  with 
restrictions  (2.7)  and  (2.9)  that  S_(s)  and  14(s)  vanish  as  s  ->  0  in  their  respective  planes 
of  analyticity.  Thus  D{s),  D+{s)  and  D-{s)  in  (2.47),  (2.56)  and  (2.57)  are  bounded  in  their 
respective  planes  of  analyticity  and  vanish  at  infinity.  Furthermore,  cr^  {x,  p)  is  expected  to 
be  square  root  singular  as  a:  — >  O'*'  and  v^{x,p)  is  expected  to  vanish  as  x  0“  to  ensure 
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continuity  of  displacement.  Prom  the  Abel  theorem  [18]  relating  asymptotic  properties  of 
transforms 


x->0+  ' 

lim  |a;|“V  (xjp)  ~ 

a;->0- 


~  lim 

S-¥<X>  '  ' 

lim  |s|^'*'’V’_(s) 

S->-00 


for  some  g  >  0.  Therefore,  it  is  expected  that  Il+{s)  ~  and  VL(s)  ~  s~^~^  as  |s|  oo, 
thus  the  products  S^.(s)L+(s)  and  L_(s)y_(s)  vanish  at  infinity.  Therefore,  each  side  of 
(2.48)  vanishes  as  |s|  -4  oo  in  the  corresponding  half  planes.  According  the  the  Liouville’s 
theorem,  a  bounded  entire  function  is  constant.  In  this  case,  W (s)  is  bounded  in  the  finite 
plane  and  W(s)  — >■  0  as  |s|  oo  so  that  the  constant  must  be  zero;  thus,  IT’(s)  =  0.  By 
using  (2.48)  and  (2.56),  the  functions  of  interest  are  then  given  by 


E+(s) 


V-{s) 


DM 

LAs) 

D-js) 

LAs) 


+  ip/Cd 

y/s-  ip/Cd  D-{s) 
P^nAs  -  w/cr)  F_(s)  ‘ 


2.2.4  Semi-infinite  smooth,  rigid  punch 


(2.58) 

(2.59) 


Consider  the  special  case  of  the  standard  semi-infinite  punch  (with  boundary  conditions 
given  by  (2.6)),  thus  a-{x,t)  =  0  and  equation  (2.42)  gives  II_(s)  =  0.  Making  further 
simplifications,  assume  the  problem  of  a  semi-infinite  smooth,  rigid  punch  such  that  fo{x,  t)  = 
fo{t).  Hence  equation  (2.43)  gives  V+(s)  =  fo{p)/is-  Therefore  (2.47)  reduces  to  D{s)  = 
LAs)fAp)/^s  and  its  factorization  follows  easily 


D{s) 

DAs) 


= m 


IS 


L-W-L-(O)  Z,-(0) 


IS 


IS 


fo(p) 


MO) 

is 


(2.60) 

(2.61) 


To  find  the  stress  intensity  factor,  an  asymptotic  expression  for  the  normal  stress  near 
the  punch  edge  is  sought.  The  Abel’s  theorem  which  relates  asymptotic  expressions  between 
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a  function  and  its  Fourier  transform  is  used  in  combination  with  the  definition  of  the  stress 
intensity  factor,  that  is,  in  the  Laplace  transform  domain 


Kjip)  —  lim  \/2'nx  a*Ax,p)  =  lim  e  S+(s). 

1-^0+  ^  s->+oo 

Substituting  (2.54)  in  (2.61)  and  the  result  in  (2.58)  it  is  found  that 

^S  +  ip/Cd 

taking  the  limit  when  s  oo  and  considering  that 

t.(0)  = 


(2.62) 


(2.63) 


(2.64) 


E+(s)  becomes 

E+(s)  =  ^^^!;(p)F40)  as  s  oo, 

(2.65) 

hence 

k;(p)  =  V2^^52^^Vp/„*(p)F-(0), 

Cr 

(2.66) 

where  .F-(O)  =  cr^ 

I  y/c^ 

1  CsCdV 

Consider  the  case  where  the  punch  approaches  the  orthotropic  half-plane  with  constant 
velocity  vq,  then  /o(p)  =  vo/p^  =  oq/ p^/^CsP^  where  ctq  —  py/^CgVQ  and  ^Jc^Cs  is  the 
dilatational  wave  speed  when  wave  propagation  is  parallel  to  the  y-axis.  Therefore 


and  by  Laplace  transform  inversion  the  dynamic  stress  intensity  factor  in  the  time  domain 
for  the  semi-infinite  punch  problem  may  be  found.  This  is 

A-,(()  =  2a„,  v/i.  (2.68) 

Note  that  (2.68)  is  the  same  as  that  for  the  case  of  the  semi-infinite  crack  problem  subject 
to  uniform  impact  loading  on  the  crack  faces  in  mode  I,  that  result  has  been  reported  by 
Rubio-Gonzalez  and  Mason  [12]. 


2.3  Finite  Punch 

2.3.1  Method  of  solution 

For  a  finite- width  punch  (0  <  a;  <  2/)  a  solution  can  be  obtained  by  superposition  of  solutions 
for  semi-infinite  punches.  Due  to  the  fact  that  any  disturbance  propagates  at  a  finite  speed, 
the  solution  for  a  semi-infinite  punch  is  valid  as  long  as  diffraction  waves  from  an  opposite 
edge  have  not  arrived  at  the  region  under  consideration. 

Consider  a  rigid  punch  of  length  21  acting  on  a  half  plane  y  >0  with  u(x,0,t)  =  fo{t) 
(figure  2.3),  the  boundary  y  =  0  outside  the  punch  is  stress  free.  The  punch  is  symmetrical, 
therefore  the  solution  to  be  found  is  symmetrical  with  respect  to  the  straight  line  x  =  I,  i.e. 

v{2l- x,y,t)  =  v{x,y,t),  (Ty{2l  -  x,y,t)  =  (Ty{x,y,t).  (2.69) 

Hence  it  is  sufficient  to  find  the  solution  in  the  region  {x  <  I,  y  >  0}. 

First  consider  the  problem  for  0  <  t  <  //c^.  Throughout  this  interval  the  diffraction 
disturbances  Mq  and  No  from  opposite  edges  of  the  punch  (re  =  0,  x  =  21)  do  not  interact, 
see  figure  2.3(a),  so  the  solution  at  x  <  I  will  be  the  same  as  in  the  problem  on  penetration 
by  a  semi-infinite  punch  (0  <  a:  <  oo).  Then  the  stress  under  the  punch  (0  <  a;  <  21)  for 
0  <t  <  l/cd  may  be  written  as 

a„(x,0,t)  =  -pv/^c,/i(«)  +  <>(*.0.<)  +  <'(2i-l,0,i),  (2.70) 

where  0,  ()  is  the  solution  for  a  semi-infinite  punch  (0  <  x  <  co)  and  -p^/^Cs/o(t) 

being  the  solution  to  the  one-dimensional  problem  on  indentation  by  an  infinite  punch 
(-00  <  X  <  oo).  Recall  that  is  the  dilatational  wave  speed  for  wave  propagation 

parallel  to  the  y-axis. 

For  Ifcd  <  t  <  2Z/cd,  figure  2.3(b),  it  follows  from  the  linearity  of  the  problem  that 
the  stress  under  the  punch  is  a  result  of  superimposing  the  disturbances  (x,  0,  t)  and 
crW(2Z  -  x,0,t)  propagating  from  the  left  and  right  edges,  respectively,  so,  equation  (2.70) 
remains  valid  throughout  this  time  interval. 
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Figure  2.3:  Schematic  representation  of  the  propagation  of  disturbances  for  (a)  0  <  t  <  l/cd-, 
(b)  l/cd  <t  <  2l/cd,  (c)  2l/cd  <t<  Sl/cd- 

For  2l/cd  <t<  Zl/cd,  a  new  pair  of  diffraction  waves  Mi  and  Ni  occurs,  figure  2.3(c). 
These  new  disturbances  must  behave  outside  the  punch  so  as  to  nullify  the  disturbances 
Mo  and  Nq  which  come  from  opposite  edges  and  violate  the  boundary  conditions  on  ?/  =  0 
outside  the  punch. 

For  2ljcd  <  t  <  31 /cd  when  these  disturbances  do  not  interact,  supplementary  distur¬ 
bances  and  within  the  domain  x  <  I  (which  are  caused  by  the  appearance  of  the 
diffraction  wave  Mi  near  the  left  edge  at  ti  >0,ti=  t-2l/cd)  are  governed  by  the  following 
boundary  conditions  on  the  axis  y  =  0: 

a«(a^>0,ii)  =  -af\2l-xM)  ^<0,  (2.71) 

0,ti)  =  0  a:  >  0.  (2.72) 

At  31 /cd  <t<  Al/cd,  disturbances  in  the  stress  under  the  punch  are  a  result  of  superpo¬ 
sition  of  the  disturbances  crW(a:,  0,ti)  and  (t^^^21  -  :r,0,ii)  caused  by  the  diffraction  waves 
Ml  and  A^i  symmetrical  with  respect  to  the  line  x  =  1.  The  same  is  valid  over  a  wider  range 
2f/cd  <t<  4l/cd. 
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Thus  once  t  exceeds  2ln/cd  (n  =  1,2,3,.. .),  there  appears  a  new  (nth)  pair  of  diffraction 
waves,  M„  and  N„.  They  carry  supplementary  disturbances  symmetrical  with  respect  to  the 
line  X  =  I  with  boundary  conditions 

=  -al^-^\2l-x,0,tn)  x<0,  (2.73) 

v^'^\x,0,tn)  =  0  x>0,  (2.74) 

with  zero  initial  conditions  at  =  0,  where  tn  =  t  —  2nl/ Cd  >  0. 

Solution  to  the  problem  with  boundary  conditions  (2.73)  and  (2.74)  is  the  same  as  in  the 
case  of  a  semi-infinite  punch  with  boundary  conditions  (2.8)  after  changing  t,  a-{x,t),  fo{x,t),  ay{x,0,t) 
andu(x,0,t)  by  -<t("-^>(2/ -  x,0,t„),  0,  cr('‘)(x,0,t„),  and  u(")(x,0,i„)  respectively.  As 
a  result,  the  stress  under  the  punch  0  <  x  <  2i  at  2kl/cd  <t<  ^l{k  +  l)/cd  {k  =  0, 1, 2, . . .) 
is 

ay{x,0,t)  =  -py/^c,m  +  E  [<)(x,0,t„)  +  c7(”)(2/  -x,0,i„)] ,  (2-75) 

n=0 

where  tn  =  t  —  2nl / Cd- 


2.3.2  Stress  intensity  factor 


An  expression  for  the  dynamic  stress  intensity  factor  near  the  punch  corners  valid  for  0  < 
t  <  Al/cd  will  be  developed  in  this  subsection.  Prom  equation  (2.75)  with  A:  =  1  we  may 


write 


0  <  t  <  2l/cd 

nJ  \  +  for  2llcd  <  t  <  4ljcd 


(2.76) 


where  Kf’\t)  and  K^^\t)  may  be  called  contributions  of  zero  and  first  orders  respectively. 
Kf\t)  corresponds  to  the  stress  intensity  factor  for  the  semi-infinite  punch,  equation  (2.68), 
this  is  associated  with  a^\x,0,t).  That  is 


(2.77) 


On  the  other  hand,  takes  into  account  the  effect  of  the  appearance  of  the  first 

diffraction  wave  Mi,  so  is  associated  with  <7(^)(x,0,t). 
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According  to  (2.71)  and  (2.72) 

E^l\s)  =  (2.78) 

Vl^\s)  =  0.  (2.79) 


.  cooiwii 


To  find  the  stress  intensity  factor  an  asymptotic 

fQ].  thg  semi-infinite  punch  is  given  by  (2.63),  then  using 


is  required.  First  note  that  E+^(s) 

F+{—s)  =  F_(s),  equation  (2.78)  becomes 

eL'^(s)  =  -e2*'*tiilMl,_(0)/o*(p)^^y^.  (2.80) 

yj-s  +  ip/cd 

From  (2.56),  (2.55)  and  (2.47)  with  I4(s)  =  0  and  E_(s)  given  by  (2.80)  it  possible  to 
show  that 

L-(0)/o*(p)  1  /•~+*^o  {-z  +  ip/cR)^z  +  iplcd  F.{z)  ^2,;,  dz  ^  ^2.81) 

27ri  J-oo+iXo  (z  +  ip/cR)^—z  -I-  ip/cd  ° 


£>+(s,p)  = 


z  —  s 


inserting  this  expression  into  (2.58),  the  desired  asymptotic  expression  is  obtained  as 

1  roo+iXo  {-z  -I-  ip/cR)^Jz  +  ip /Cd 


y/—ty/s  Cr  (^z  +  Ip/cr)‘\J—Z ip/Cd 


F-lz)  Ozlij 

-^r-^^e  dz  as  s  -»  oo 
zF^iz) 


(2.82) 


where  equation  (2.64)  for  L_(0)  has  been  used.  Making  the  change  of  variable  ^  =  ipC  and 
using  (2.62),  the  stress  intensity  factor  in  the  Laplace  transform  domain  may  be  written  as 


K 


;“>(?)  =  -V2iin(^VPK(p)f'-(D)^  r*’"" me-^'di;  (2.83) 

Cr  /TTl  Jtq—ioo 


where 


(-C  +  l/cit)v/C  + 1/crf  F.{0 

(C  +  i/cr)\/~C  +  1/^d  C^+iO 


(2.84) 


To  perform  the  inversion  of  Kj{p)  the  Cagniard-de  Hoop  method  can  be  applied.  The 
central  idea  of  the  Cagniard-de  Hoop  scheme  [20,  16]  is  to  convert  the  integral  in  (2.83)  to  a 
form  which  allow  inversion  of  the  one-sided  Laplace  transform  by  observation.  The  path  of 
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integration  is  modified  to  form  a  closed  contour,  figure  2.4,  such  that  the  integrand  in  (2.83) 
is  analytic  inside  of  this  contour.  Consider  again  the  case  where  the  punch  approaches 
the  orthotropic  half-plane  with  constant  velocity  Vq,  then  /o(p)  =  Vo/p^  =  (^ol Py/^Cgp^. 
Applying  Cauchy’s  theorem,  Jordan’s  lemma  and  the  fact  that  /(C)  =  /(C)  can  write 
as 


AT, •<')(?)  =  r  Iml/(C)|e-«'<iC. 

I  yy/  A-  s  py/c^Cs-K  '  p^'^  Jl/cd 


(2.85) 


Letting  2C/  =  this  equation  becomes 

rim|/(>l/201e-’"’/f(>)  -  il/cMn-  (2.86) 

^  ^  Cr  py/^Csn2l  p^/^  Jo 

and  the  inversion  of  the  Laplace  transform  becomes  obvious.  Equation  (2.86)  is  a  product 
of  two  transforms,  so  that  Kj^\t)  is  a  convolution  of  the  inverse  of  the  two  transforms,  i.e. 

*:}'>(()  = /'  dn,  (2.87) 

CRy/^  TT  V  TT  Jl/c^ 


and,  by  substitution  of  F-{0)  yields 


ArP>(t)  = 


^Cs  f* 

\  y/^  Jl/Cd 


f  lm[I(‘n)]^Ji  -  T]  dr]. 

j\Ica 


where  i  =  1/21 . 


(2.88) 


l/cd  1/cr\  ReQ 


Figure  2.4:  Contour  of  integration  to  evaluate  the  integral  in  (2.83)  in  the  C—  plane. 
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2.4  Semi-infinite  Crack 


The  semi-infinite  crack  problem  is  considered  in  this  section.  The  analysis  follows  that  of 
the  punch  case  with  some  minor  differences.  So  there  is  no  need  to  repeat  the  analysis  of 
the  last  sections  and  only  the  important  points  will  be  outlined. 

2.4.1  Method  of  Solution 

Consider  the  plane  problem  of  an  infinite  orthotropic  medium  containing  a  semi-infinite 
crack,  figure  2.5.  The  crack  faces  are  along  the  negative  a;— axis  and  the  origin  of  the  xy 
coordinate  system  is  in  the  crack  tip.  Let  a_(a:,  t)  be  the  pressure  applied  on  the  crack  faces, 
and  consider  a  more  general  problem  where  the  displacement  ahead  of  the  crack  tip  is  given 
by  fo{x,t),  thus  the  boundary  conditions  are 

ay{x,0,t)  =  a-.{x,t)  for  —  oo  <  a;  <  0, 

Txy{x,0,t)  =  0  for  —  oo  <  a;  <  CO,  (2.89) 

v{x,0,t)  =  for  a:  >  0. 

In  addition,  the  condition  of  zero  displacements  at  infinity  and  zero  initial  conditions  are 
assumed.  Note  that  with  fo{x,t)  =  0  the  standard  problem  with  the  semi-infinite  crack  is 
recovered  [12].  The  function  /o(a;,  t)  is  assumed  to  be  zero  at  any  x  >  Cd,t  +  k,  (k  >0). 

o„H(t)  . . . 

OcHd) 

Figure  2.5:  Schematic  of  the  semi-infinite  crack  geometry. 

The  differential  equations  describing  this  problem  are  (2.1)  and  (2.2)  with  boundary 
conditions  (2.89).  Therefore  the  same  procedure  as  for  the  punch  problem  is  applied,  i.e.. 
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application  of  Laplace  and  Fourier  transforms  to  the  equations  of  motion  and  boundary 
conditions  following  by  the  use  of  the  Wiener-Hopf  technique  to  find  ^+(5)  and  VL(5)  which 
are  given  by  (2.58)  and  (2.59).  It  should  be  noted  that  in  this  case,  the  planes  of  analyticity 
of  the  double  transforms  11+ (s)  and  V_(s)  will  be  Im(s)  >  —p/cd  and  Im(s)  <  0  respectively. 

For  the  standard  semi-infinite  crack  problem  /o(iC,  t)  =  0  and  cr_(x,  t)  =  —aoH[t),  hence 
y+(s)  =  0  and  E_(s)  =  -ao/ips.  Therefore  (2.47)  reduces  to  D{s)  =  -L+(s)cro/zps  and  its 
factorization  follows  easily 


D{s) 

D-{s) 


L+(s)  -  L+(0)  ^  L+(0) 

is  is 

<7o  L+(0) 

p  is  ’ 


ips  p 


(2.90) 

(2.91) 


and  hence  the  unknown  functions  11+ (s)  and  ^-(s)  are 


E+(s) 

V.{s) 


{s  +  ip/cR)ao 
ijs  +  ip/cd  P 
y's  -  ip/cR 
-  ip/cd) 


L+{s)  —  ■L+(0) 
is 

ap  L+  (0)  1 

P  is  F+(s)' 


his), 


(2.92) 

(2.93) 


Following  the  same  procedure  as  in  the  semi-infinite,  smooth,  rigid  punch  problem  we 
can  get  the  stress  intensity  factor  for  the  standard  semi-infinite  crack  problem  with  uniform 
impact  pressure  on  the  crack  faces,  this  is 


Kj(t)  =  2<7o^ 


Vt. 


(2.94) 


2.5  Finite  Crack 


2.5.1  Method  of  Solution 


A  solution  for  the  finite  crack  problem  will  be  obtained  by  superposition  of  solutions  for 
semi-infinite  cracks.  Consider  an  infinite  orthotropic  plate  containing  a  crack  of  length  21 
with  crack  faces  parallel  to  the  re— axis  and  the  crack  tips  located  at  rc  =  —21  and  re  =  0, 
figure  2.6.  By  symmetry  it  is  expected  that  the  solution  is  symmetric  with  respect  to  the 
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line  X  =  —I,  that  is 


v{-2l-x,y,t)  =  v{x,y,t),  ay{-2l  -  x,y,t)  =  (Jy{x,y,t).  (2.95) 

Hence  it  is  sufficient  to  find  the  solution  in  the  region  {x>  —I,  y  >  0}. 

For  the  time  interval  0  <  t  <  2l/cdi  the  stress  field  is 

ay{x,  0,  t)  =  -ao  +  af{x,  0,  t)  +  af  (-2/  -  x,  0,  t),  (2.96) 

where  (7^°^(x,0,t)  is  the  solution  for  a  semi-infinite  crack. 

For  21 /cd  <  t  <  Al/cd,  supplementary  disturbances  and  within  the  domain 
X  >  —I  (which  are  caused  by  the  appearance  of  the  diffraction  wave  Ni  near  the  right  tip  at 
ti>  0,ti  =  t-  21 /cd)  are  governed  by  the  following  boundary  conditions  on  the  axis  y  =  0: 

cr^^^(a;,0,ti)  =  0  x<0,  (2.97) 

0,  ti)  =  —v^^\—2l  —  x,0,ti)  x  >  0.  (2.98) 

Note  that  the  boundary  condition  (2.98)  is  just  to  nullify  the  disturbances  that  come  from 
opposite  crack  tips  and  violate  the  boundary  conditions  on  y  =  0  outside  the  crack. 


Figure  2.6:  Schematic  representation  of  the  propagation  of  disturbances  for  (a)  0  <  t  <  Z/cd, 
(b)  l/cd<t  <  21 /cd,  (c)  21 /cd  <t  <31/ Cd- 

Thus  once  t  exceeds  2ln/cd  {n  =  1, 2, 3, . . .),  there  appears  a  new  (nth)  pair  of  diffraction 
waves,  Mn  and  Nn-  They  carry  supplementary  disturbances  symmetrical  with  respect  to  the 
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line  X  =  —I  with  boundary  conditions 


(2.99) 

(2.100) 


crJ”H^.0.*n)  =  0  a;<0, 

u(")(a:,0,<„)  =  -u("-i)(-2Z-x,0,<n)  a:  >  0, 


with  zero  initial  conditions  at  =  0,  where  t  —  2nl/cd  >  0. 

Solution  to  the  problem  with  boundary  conditions  (2.99)  and  (2.100)  is  the  same  as  in  the 

case  of  a  semi-infinite  crack  with  boundary  conditions  (2.89)  after  changing  t,  a-{x,t),  /o(a:,  t),  cFy{x,0,t^ 

and  7;(a:,0,t)  by  0,  -u(”-i>(-2/  -  a;,0,t„),  a(”^x,0,t„),  and  u(”)(a;,0,t„)  respectively. 

As  a  result,  the  stress  field  at  2kl/cd  <t<  2l{k  -f  l)/crf  {k  =  0, 1, 2, . . .)  is 

ay{x,  0,  t)  =  -<7o  -t-  E  0,  tn)  +  aW(-21  -  rc,  0,  Q]  ,  (2.101) 

n=0 

where  tn  =  t  —  2nl/cd. 


2.5.2  Stress  intensity  factor 

An  expression  for  the  dynamic  stress  intensity  factor  near  crack  tips  valid  for  0  <  t  <  4// Cd 
is  developed  next.  From  equation  (2.101)  with  A:  =  1  we  may  write 

X  (f\  —  {  0  <  t  <  21 /cd  /2  ^Q2’j 

for  2l/cd<t<4lfcd 

being  Kf\t)  and  the  contributions  of  zero  and  first  orders  respectively. 

corresponds  to  the  stress  intensity  factor  for  the  semi-infinite  crack,  equation  (2.94),  this  is 
associated  with  (r(°)(rc,0,t).  That  is 


(2.102) 


Af  (t)  =  2(70, 


TTy^^ 


(2.103) 


On  the  other  hand,  Kj^\t)  takes  into  account  the  effect  of  the  appearance  of  the  first 
diffraction  wave  Ni,  so  Ki^\t)  is  associated  with  cr^^)(a:,  0,  t). 

According  to  (2.97)  and  (2.98) 

=  0, 


vr\s) 


-e-2*'Vio)(-s). 
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(2.104) 

(2.105) 


Figure  2.7:  Stress  intensity  factor  history  for  a  finite  crack  in  boron-epoxy  composite  type  1. 

To  find  the  stress  intensity  factor  an  asymptotic  expression  of  as  s  ^  oo 

is  required.  First  note  that  for  the  semi-infinite  crack  is  given  by  (2.93),  then  using 

F+{-s)  =  F_{s),  equation  (2.104)  becomes 

_ ^MO).  (2.106) 

^^\2^{-s-ip|cd)pF+{s)  is 

From  (2.56),  (2.55)  and  (2.47)  with  E_(s)  =  0  and  V+(s)  given  by  (2.106)  it  possible  to 
show  that 


D+is,p)  = 


<toL+(0)  1  {z-ip/cR)^-z-ip/cd  F-{z) 

•  2m  J-  '  -  ■  '  ‘ 


oo+tAo  (—2;  —  iplcR)\jz  —  ip/cd 


z  —  s' 


(2.107) 


inserting  this  expression  into  (2.58),  the  desired  asymptotic  expression  is  obtained  as 


Wf  -  1  1  (^  -  iplCR)y/-Z  -  ip/Cd 

^  iy2y^F+(0)  y/sp^^^  2m  7-00+1  Ao  (—z  —  iplcR)sJz  —  ip/cd 


F-  {z)  _2zii  , 

—r-^e  dz  as  s  ->  00 
zF^{z) 


(2.108) 


where  T+(0)  =  cr/ y/^^dF+{Qi)  has  been  used.  Making  the  change  of  variable  2:  =  -ipC 
and  using  (2.62),  the  stress  intensity  factor  in  the  Laplace  transform  domain  may  be  written 
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Finite  Crack,  Boron-epoxy  type  I 


— 1 - 1 - 1 - 1 

6X3CI  SOI* 

-  using  FFT 

_ I _ I _ ^ ^ - L_J - 1 - ^ - . - 1 - 1 - 1 

0  2  4  6  8  10  12  14  16  18  20 

Cdt/2I 


Figure  2.8:  Comparison  between  the  exact  solution  of  Ki{t)  and  the  results  derived  using  the 
method  of  Kassir  (1983)  but  with  FFT  for  the  Laplace  transform  inversion,  for  boron-epoxy 
composite  type  I. 


as 


where 


V^apCR  1 _ 1_ 

27ri 


rro+ioo 
J  TO  —ZOO 


(C  +  l/cji)y-C+ 1/cd  F+{C) 

(— c  +  ^/cr)\Jc + 1/crf  c-^-(c) 


(2.109) 


(2.110) 


To  perform  the  inversion  of  Kj{p)  the  Cagniard-de  Hoop  method  can  be  applied  as  in 
the  punch  problem.  After  doing  the  inversion  we  get  the  first  order  contribution  to  the 
dynamic  stress  intensity  factor  in  the  time  domain,  as 


Kf\t)  = 


TT^ 


[  lm[I{ri)]\li-r]dr]. 
y'C22  Jl/cd 


(2.111) 


where  i  =  t/2l  and  /(C)  is  given  by  (2.110). 
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K(t)  Finite  Punch,  Graphite-Epoxy  Composite 


K(t)  Finite  Punch,  Isotropic  Material 


0  2  4  6  8  10  12 

c.t/21 


(c) 


Figure  2.9:  Dynamic  stress  intensity  factor  for  the  punch  problem,  (a)  graphite-epoxy  com¬ 
posite,  (b)  graphite-epoxy  composite  with  fibers  parallel  to  the  y-axis,  (c)  isotropic  material. 


2.6  Results  and  Conclusions 

2.6.1  Punch  Problem 

Figure  2.9  illustrates  the  dynamic  stress  intensity  factor  history  Ki{t)  for  the  punch  problem 
for  different  materials  listed  in  Table  2.1.  Note  that  Ki{t)  for  the  finite  punch,  equation 
(2.76),  is  valid  only  in  0  <  Cdt/2f  <  2,  however  it  has  been  plotted  for  an  extended  time 
interval  to  show  that  the  first  order  contribution  K\^\t)  is  small  even  outside  the  interval 
of  validity  0  <  Cdt/2l  <  2.  Figure  2.9(b)  shows  that  for  the  graphite-epoxy  composite  with 
fibers  parallel  to  the  y-axis  the  stress  intensity  factor  departs  from  the  semi-infinite  punch 
result  at  earlier  normalized  time  than  in  figure  2.9(a)  where  the  fibers  are  along  the  x— axis 
direction.  The  stress  intensity  factor  is  normalized  with  respect  to  Kq  =  aoVTd- 


Finite  Crack.  Boron-epoxy  type  li 


Figure  2.10:  Comparison  between  the  exact  solution  of  Ki{t)  and  the  results  derived  using 
the  method  of  Kassir  (1983)  but  with  FFT  for  the  Laplace  transform  inversion,  for  boron- 
epoxy  composite  type  II. 

2.6.2  Crack  Problem 

Figures  2.7-2.11  show  the  dynamic  stress  intensity  factor  Kj{t)  for  the  finite  crack  in  different 
materials.  The  labels  P  and  R  represent  the  arrival  time  of  the  dilatational  and  Rayleigh 
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waves  respectively  from  the  opposite  crack  tip.  The  label  2P  means  the  time  needed  for 
the  dilatational  wave  to  travel  the  crack  length  twice,  i.e.,  the  solution  derived  here  for  the 
finite  crack  is  valid  until  the  time  labeled  with  2P.  The  curved  marked  with  ’’using  FFT” 
in  figures  2.7  and  2.8  represent  the  result  obtained  for  Kj{t)  using  the  method  of  Kassir 
and  Banyopadhyay  [13] ,  but  with  different  numerical  schemes  to  solve  the  Fredholm  integral 
equation  and  using  fast  Fourier  transform  for  Laplace  transform  inversion  [21].  Comparing 
these  results  with  the  exact  solution  just  developed  it  is  concluded  that  such  schemes  were 
more  appropriate  as  was  claimed  in  [21].  The  dynamic  overshoot  is  better  described.  In 
all  the  cases,  figures  2.7  and  2.8,  the  results  obtained  by  Kassir  and  Banyopadhyay  [13]  are 
shown  to  depart  from  the  exact  solution.  Figure  2.11  shows  Kiit)  for  finite  crack  in  graphite- 
epoxy  composite  with  fibers  parallel  to  the  y— axis,  in  this  case  the  dynamic  overshoot  occurs 
inside  the  time  interval  of  validity  of  the  exact  solution,  i.e.,  in  0  <  Cdt/2l  <  2,  and  hence  it 
is  accurately  described.  The  normalization  factor  is  the  long  time  stress  intensity  factor  for 
a  finite  crack,  Ki(oo)  =  <to\/^. 


Figure  2.11:  Stress  intensity  factor  history  for  a  finite  crack  in  graphite-epoxy  composite 
with  fibers  parallel  to  the  y— axis. 
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Graphite 

Epoxy 

E-Glass 

Epoxy 

Boron 

Epoxy 

Epoxy 

Cll 

20.77 

C22 

2.18 

3.12 

C12 

0.49 

0.52 

0.79 

A*i2  (GPa) 

7.48 

5.5 

6.4 

p  (Kg/m^) 

1600 

2100 

1990 

1260 

Table  2.1:  Mechanical  properties  used  for  the  analysis. 
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Chapter  3 

Experimental  Investigation  of 
Dynamic  Punch  Tests  on  Isotropic 
and  Composite  Materials 
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Overview 


An  experimental  investigation  is  conducted  on  the  two-dimensional  punch  problem  for  isotropic 
materials  and  unidirectional  fiber-reinforced  composite  materials  under  quasi-static  and  im¬ 
pact  loading.  Singular  stresses  are  generated  in  the  specimen  near  the  punch  corners  and  the 
stress  intensity  factor  Kj  is  introduced  to  describe  the  singular  stress  field.  Laser  interferom¬ 
etry  is  used  to  measure  in-plane  stresses  (transmission  mode)  and  out-of-plane  displacements 
(reflection  mode),  and  then  estimate  the  stress  intensity  factor.  In  the  dynamic  case,  a  high 
speed  photography  technique  was  employed  to  capture  the  transient  response  of  the  speci¬ 
men  and  measure  K{t)  just  after  the  impact.  In  all  the  cases  a  good  agreement  between  the 
measurements  of  K  and  the  theoretical  predictions  was  found. 

3.1  Introduction 

The  potential  use  of  composite  materials  in  the  construction  of  defense  structures  makes  it 
necessary  to  understand  the  initiation  of  damage  in  the  composite  as  caused  by  impacting 
fragments.  Typically,  this  kind  of  damage  initiation  is  simulated  in  the  laboratory  by  the 
punch  test.  In  such  a  test,  fiber-reinforced  epoxy  matrix  composites  will  develop  a  limited 
amount  of  plasticity  before  failure  occurs.  This  fact  suggests  an  elastic  solution  of  the  stress 
field  will  play  an  important  role  in  understanding  and  predicting  the  behavior  of  composites 
under  punch  test,  impact  conditions. 

Like  in  fracture  mechanics,  a  stress  intensity  factor  Kj  is  introduced  in  this  work  to 
characterize  the  behavior  of  the  singular  stresses  around  the  punch  corners.  This  parameter 
is  measured  experimentally  using  laser  interferometry  and  then  compared  with  theoretical 
predictions.  The  experimental  investigation  is  conducted  on  the  quasi-static  and  dynamic 
punch  test  for  isotropic  and  orthotropic  materials.  In  the  test,  a  rigid  punch  comes  in  contact 
with  one  side  of  the  specimen,  figure  3.1.  The  lateral  shearing  interferometer  of  coherent 
gradient  sensing  (CGS),  in  conjunction  with  high  speed  photography,  is  used  to  obtain  real 
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time  interferograms  of  the  singular  stress  field  generated  near  the  corner  of  the  punch. 


The  application  of  the  CGS  method  in  the  punch  test  is  new  and  the  verification  of 
the  dynamic  solution  is  new  as  well.  Here,  results  for  isotropic  as  well  as  for  composite 
materials  are  reported.  In  the  isotropic  case,  the  material  selected  was  PMMA  (polymethyl 
methacrylate)  a  transparent  acrylic  for  which  the  CGS  method  in  transmission  mode  is  well 
suited.  For  the  orthotropic  case,  the  material  chosen  was  a  unidirectional  graphite-epoxy 
composite  with  fiber  volume  fraction  of  0.65. 

The  CGS  method  has  been  applied  successfully  in  dynamic  fracture  experiments  on 
isotropic  materials  [1,  2]  and  more  recently  it  has  been  applied  to  dynamic  fracture  of  fiber 
reinforced  composite  materials  by  Lambros  and  Rosakis  [3,  4]  and  in  the  analysis  of  the 
quasi-static  fracture  behavior  of  composites  by  Liu  et  al.  [5].  The  advantage  of  the  CGS 
method  combined  with  high  speed  photography  is  that  real-time  visualization  of  dynamic 
fracture  initiation  and  growth  can  be  obtained  and  measurements  of  the  dynamic  stress 
intensity  factor  during  these  processes  can  be  made. 

The  CGS  method  is  applied  to  measure  the  stress  intensity  factor  on  the  punch  test  on 
four  sets  of  experiments;  two  static  and  two  dynamic.  In  both  cases  the  materials  employed 
for  the  specimens  are  PMMA  and  graphite-epoxy  composites.  In  the  quasi-static  test,  K[ 
is  compared  successfully  with  known  solutions  of  the  punch  problem  for  isotropic  and  or¬ 
thotropic  materials.  In  the  dynamic  case,  each  photograph  taken  by  the  high  speed  camera 
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( 


is  digitized  and  Ki{t)  is  then  measured.  The  time  evolution  of  Kj  agrees  very  well  with 
the  theoretical  predictions  for  the  dynamic  punch  problem  developed  in  Rubio-Gonzalez  and 
Mason  [6]. 

3.2  Elastodynamic  Analysis  of  the  Finite  Punch  Prob¬ 
lem 


The  punch  problem  is  of  great  importance  in  solid  mechanics  for  its  multiple  technical 
applications  including  ballistic  impact,  metal  forming  and  manufacturing  operations  such 
as  punching  and  blanking.  While  the  quasi-static  punch  problem  is  well  studied  in  contact 
mechanics  [7,  8],  there  is  little  work  on  the  dynamic  case.  In  the  dynamic  case  it  is  assumed 
that  the  punch  approaches  the  material  with  a  fixed  velocity,  and  that  wave  propagation  is 
involved  thereby  complicating  the  mathematical  analysis. 

As  shown  by  Roessig  and  Mason  [9]  and  by  Rubio-Gonzalez  and  Mason[6]  the  dynamic 
stress  intensity  factor  for  the  two-dimensional  dynamic  semi-infinite  punch  coincides  with 
that  for  the  semi-infinite  crack  under  uniform  impact  loads,  cxq,  applied  on  the  crack  faces 
in  mode  I.  An  expression  for  for  the  dynamic  stress  intensity  factor  has  been  developed  by 
Rubio-Gonzalez  and  Mason  [10],  it  is  given  by 

Kr{t)  =  2ao.  Vi,  (3.1) 

\  T^y/C22 

where 


and  Cs  =  represents  the  velocity  of  the  in-plane  shear  wave  propagating  along  the 

the  principal  material  axes  and  p  is  the  mass  density.  The  non-dimensional  constants  Cjj  are 
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given  by 


El 

C22  =  {E2/E\)C\i, 

C12  =  ,  1^12^22  =  ^^21^11, 

for  generalized  plane  stress,  and  by 

C12  =  - x(^21  +  -^*^13*^32)) 

Ati2A  £;i 

A  =  1  —  1^121^21  ~  t^2Z^Z2  —  ^Zl^lZ  ~  Vi2V22Vz\  —  1^13*^21*^32) 


(3.4) 


(3.5) 


for  plane  strain.  In  these  equations  Ei,  Hij  and  i/ij  {i,j  =1,2,3)  are  the  engineering  elastic 
constants  of  the  material  where  the  indices  1,  2,  and  3  correspond  to  the  directions  {x,y,  z) 
of  a  system  of  Cartesian  coordinates  chosen  to  coincide  with  the  axes  of  material  orthotropy. 

The  problem  of  a  finite,  rigid  and  flat  punch  impacting  an  orthotropic  half-plane  has 
been  analyzed  by  Rubio-Gonzalez  and  Mason  [6]  for  the  orthotropic  case  using  integral 
transforms  and  the  Wiener-Hopf  technique  and  by  Roessig  and  Mason  [9]  for  the  isotropic 
case.  A  solution  for  the  dynamic  stress  intensity  factor  Kj{t)  was  developed  which  is  valid 
while  the  dilatational  wave  travels  the  punch  width  twice.  That  is 


Ki 


for 

Kf\t)  +  Ki^\t)  for 


0  <  i  <  2l/cd 
2l/cd  <t  <  Al/cd 


(3.6) 


where  Kf\t)  and  Kj^\t)  may  be  called  contributions  of  zero  and  first  orders  respectively. 
Kf\t)  corresponds  to  the  stress  intensity  factor  for  the  semi-infinite  punch,  equation  (3.1). 
On  the  other  hand,  KY^{t)  takes  into  account  the  effect  of  the  appearance  of  the  first 
diffraction  wave  emitted  by  the  opposite  punch  corner.  The  velocity  Cd  =  y/^Cg  represents 
the  dilatational  wave  speed  along  the  x— axis. 
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It  was  found  that  Kj^\t)  makes  a  small  contribution  to  Kj{t)  in  equation  (3.6).  Thus, 
for  the  dynamic  punch  test,  the  experimental  results  will  be  compared  to  given  by 

equation  (3.1). 

It  has  been  shown  [10]  that  the  orthotropic  formulation  includes  the  isotropic  result 
as  special  case.  That  is,  Ki{t)  for  isotropic  materials  is  recovered  from  the  orthotropic 
expression  (3.1)  with  the  proper  substitution  of  the  elastic  constants  and  it  is  given  by 

Ki{t)  =  (3.7) 

{l-v) 

This  expression  was  developed  for  the  isotropic  case  by  Freund.  [11] 

Note  that  in  both  cases,  equations  (3.1)  and  (3.7),  the  stress  intensity  factor  is  propor¬ 
tional  to  the  square  root  of  time  and  may  be  written  in  the  form 

Ki{t)  =  cToC/x/t, 


where  the  material-related  constant  Ci  is 

Ci  =  2, 

for  the  orthotropic  case,  and 


2cs^ 


(3.8) 


c,  =  -  ci/cl)/ir,  (3.9) 

for  the  isotropic  case. 

Equations  (3.1)  and  (3.7)  were  derived  under  the  assumption  that  the  semi-infinite  crack 
face  loading  (or  semi-infinite  punch)  was  suddenly  applied.  Once  Ki{t)  is  established  for  the 
step  loading,  the  corresponding  result  for  any  time  loading  may  be  obtained  by  superposition 
[11].  For  the  particular  case  of  a  step  load  of  finite  duration  t*,  the  resulting  stress  intensity 
factor  is  given  by 

Ki{t)  =  aoCi  -  t*)] ,  (3.10) 

where  Ci  is  given  by  equations  (3.8)  and  (3.9)  for  the  orthotropy  and  isotropic  cases,  respec¬ 
tively. 
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3.3  Coherent  Gradient  Sensing  (CGS) 


CGS  is  a  full  field  optical  method  that  may  be  used  to  study  deformations  in  transparent  as 
well  as  opaque  solids.  This  technique  provides  measures  of  either  in-plane  gradients  of  stresses 
(transmission)  or  in-plane  gradients  of  out-of-plane  displacements  (reflection).  Consider 
a  specimen  of  uniform  thickness  in  the  undeformed  state  with  the  mid-plane  coinciding 
with  the  a:  -  y  plane,  see  figure  3.2.  The  specimen  is  loaded  so  that  it  causes  nonuniform 
spatial  gradients  in  the  optical  path  when  light  is  transmitted  through  it  or  reflected  from  its 
surface.  These  gradients  are  related  to  the  deformation.  For  transmission  of  light  through  a 
transparent  specimen  gradients  in  refractive  index  are  related  to  gradients  of  stress  through 
the  stress-optic  coefficient,  for  reflection,  gradients  in  nonuniform  out-of-plane  displacements 
are  related  to  deformation  through  the  Poisson  effect. 

Experimental  Arrangement 

The  schematic  of  the  experimental  setup  used  for  transmission  CGS  is  shown  in  figure  3.2(a). 
A  transparent  optically  isotropic  plate  specimen  is  illuminated  by  a  collimated  beam  of 
coherent  laser  light.  The  transmitted  wave  passes  through  a  pair  of  high  density  gratings 
Gi  and  G2,  separated  by  a  distance  A.  The  field  distribution  on  the  G2  plane  is  spatially 
transformed  by  the  lens  Li.  By  locating  a  filtering  aperture  around  either  the  ±1  diffraction 
orders  at  the  filter  plane,  information  regarding  the  stress  gradients  is  obtained  on  the  image 
plane  of  the  lens  L2. 

Figure  3.2(b)  shows  the  modification  of  the  above  setup  for  measuring  surface  deflection 
of  opaque  solids  when  studied  in  reflection  mode.  In  this  case,  the  specularly  reflecting 
specimen  surface  is  illuminated  by  a  collimated  beam  of  laser  light.  The  reflected  beam,  as 
in  the  previous  case,  is  processed  through  the  optical  arrangement  which  is  identical  to  the 
one  shown  in  figure  3.2(a). 

A  brief  discussion  of  the  analysis  of  the  system  as  described  by  Mason  et  al.  [12]  now 
follows.  After  exiting  from  the  deformed  specimen,  the  beam  falls  upon  the  first  of  two 
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Figure  3.2:  Schematic  of  the  experimental  setup  for  (a)  transmission  CGS,  and  (b)  reflection 
CGS. 

identical  diffraction  gratings,  see  figure  3.3.  The  primary  grating  Gi  splits  the  beam  into  a 
direct  beam  and  numerous  diffraction  orders.  For  the  sake  of  brevity  only  the  first  diffraction 
orders  (±1)  and  the  direct  beam  are  considered.  The  second  diffraction  grating  diffracts  both 
the  direct  beam  and  the  first  diffraction  orders  into  three  beams  each,  giving  a  total  of  nine 
beams  behind  the  second  grating.  Of  these  nine  beams  the  (0,±1)  and  (±1,0)  are  parallel. 
A  filtering  lens  can  then  be  used  to  focus  these  parallel  beams  to  a  single  spot.  In  the 
focal  plane,  a  spatial  Alter  is  used  to  eliminate  all  superfluous  beams.  Another  lens  is  then 
required  to  produce  the  final  image. 

The  wave  that  impinges  on  the  first  diffraction  grating  Gi  is  assumed  to  be  parallel  with 
some  phase  difference,  S{x,y),  due  to  the  deformations  in  the  specimen.  The  two  gratings 
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Grating  Gl  Grating  G2  Filtering  Lens  Ll  Filter  Plane 

Figure  3.3:  Schematic  of  the  principle  of  CGS. 


shift  one  beam  with  respect  to  the  other  by  a  distance 

€  =  Atan0«A0,  (3-11) 

where  A  is  the  separation  between  the  gratings  and  6  is  the  angle  of  diffraction  as  shown  in 
figure  3.3.  The  diffraction  angle  is  given  by  the  relation 

6  =  sin”’^  —  fa  — ,  (3-12) 

P  P 

where  A  is  the  wavelength  of  the  light  and  p  is  the  grating  pitch. 

The  two  parallel,  sheared  wavefronts  constructively  interfere  if  their  phase  difference  is 
an  integer  multiple  of  the  wavelength,  or 

S{x  +  e,y)-S  (x,  y)  =  mX.  (3.13) 


Dividing  by  e  yields 

S{x  +  e,y)-  S{x,  y)  _  mX 
-  ^ 

€  e 

which,  for  small  e,  can  be  approximated  as 

dx  A 

Equation  (3.15)  is  valid  for  the  derivative  in  either  the  x  or  y  directions. 


(3.14) 


(3.15) 
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For  a  transparent  materials,  the  phase  difference  S{x,y)  is  composed  by  two  contribu¬ 
tions;  the  changes  in  specimen  thickness  due  to  Poisson’s  contraction,  and  the  change  in 
refractive  index  of  the  material  due  to  variations  in  hydrostatic  stress  (Maxwell  relation).  It 
is  possible  to  show  [1]  that  for  an  isotropic  and  linearly  elastic  material  under  plane  stress 
conditions,  S{x,y)  is  given  by 

S{x,  y)  «  ch{(Jx  +  cFy),  (3-16) 


where  c  is  a  modified  stress-optic  constant,  h  is  the  thickness  and  cTx  S'Hd  (^y  plane  stress 
thickness  averages  of  the  stress  components  in  the  material.  Substituting  (3.16)  into  (3.15) 
gives  the  result 


<*£(2^  =  !^,  „.  =  o,±i,±2.... 

ax  A 


(3.17) 


For  an  opaque  material  reflecting  the  incident  laser  light,  the  phase  difference,  S{x,y), 
is  wholly  attributed  to  changes  in  specimen  thickness  due  to  lateral  contraction,  thus 


Six,y)  =  2w,  (3.18) 

where  w  is  the  z— component  of  the  displacement  vector.  Substituting  (3.18)  into  (3.15) 
yields 


dw 

dx 


mp 


m  =  0,  ±1,  ±2, . . .  (3.19) 

Summarizing,  the  formation  of  constructive  interference  on  the  image  plane  is  governed 


by 

transmission 


reflection 


d{ax  +  (Ty) 
Cfl  ^ 

ox 

= 

mp 

X’ 

m  =  0,  ±1,  ±2, . . . 

(3.20) 

j^d{ax  +  cry) 

ch 

oy 

= 

np 

n  =  0,  ±1,  ±2, . . . 

(3.21) 

dw 

dx 

mp 

m 

=  0,±1,±2,... 

(3.22) 

dw 

dx 

np 

n  ’ 

=  0,±1,±2,... 

(3.23) 
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where  h  is  the  nominal  thickness  of  the  specimen,  p  the  grating  pitch,  A  the  separation 
between  gratings  and  c  the  modified  stress-optic  constant. 

3.4  Preparation  and  Characterization  of  the  Compos¬ 
ite  Specimens 

The  preparation  and  characterization  of  the  specimens  of  PMMA  demanded  a  minimum 
of  work;  however  the  preparation  and  characterization  of  composite  specimens  required  a 
careful  procedure  in  order  to  successfully  apply  the  CGS  technique.  Details  for  the  latter 
case  are  given  below. 

3.4.1  Tension  Test  of  Composites 

The  graphite-epoxy  composite  plates  used  in  the  experiments  were  provided  by  Composite 
Mirror  Applications  Inc.  of  Tucson,  Arizona. 

First,  a  series  of  tension  tests  were  performed  to  determine  the  material  properties  of 
this  particular  graphite-epoxy  composite.  The  material  was  assumed  to  be  transversely 
isotropic.  The  following  properties  were  measured  :  Young’s  modulus  Ei,  E2,  Poisson’s 
ratio  V12  and  shear  modulus  /X12  (see  figure  3.4  for  definition  of  axes).  The  tension  test  was 
conducted  according  with  recommendations  of  the  ASTM  standard  D  3039M-95a  Standard 
Test  Method  for  Tensile  Properties  of  Polymer  Matrix  Composite  Materials  and  performed 
on  screw  driven  machine. 

Table  3.2  summarizes  the  material  parameters  for  graphite-epoxy  composite  used  in  the 
experiments.  In  table  3.2  some  calculated  material  parameters  are  given. 

3.4.2  Specimen  Preparation 

As  mentioned  earlier,  an  optically  flat  and  specularly  reflecting  surface  is  needed  in  order  to 
apply  the  CGS  technique  in  reflection  mode.  Therefore  a  very  high  quality  surface  prepa¬ 
ration  is  required  on  the  composite  material.  Due  to  the  microstructure  of  the  composite, 
it  is  not  possible  to  directly  deposit  a  reflective  film  on  the  specimen  by  vacuum  deposition. 
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Table  3.1:  Material  properties  of  the  particular  graphite-epoxy  composite  used  in  the  exper- 


El  (GPa) 

92  ±3 

(GPa) 

6.7  ±0.1 

1^12 

0.38  ±0.01 

fii2  (GPa) 

3.00  ±0.05 

P  (Kg/m^*) 

1460  ±20 

Table  3.2:  Calculated  material  parameters  of  the  particular  graphite-epoxy  composite  used 
in  the  experiments.  Plane  stress  conditions _ _ 


Cll 

30.85 

Cl2 

0.85 

C22 

2.25 

^12  (GPa) 

3.00  ±0.05 

p  (Kg/m^) 

1460  ±20 

Cs  (m/s) 

1430  ±12 

Cd,x  (m/s) 
parallel  to  the  fibers 

7970  ±13 

Cd,y  (m/s) 

perpendicular  to  the  fibers 

2160  ±16 

1/23  (m/s) 

.30  ±0.01 
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Figure  3.4:  Micrographs  of  the  unidirectional  graphite-epoxy  composite  with  fibers  along 
the  a;— axis. 

A  procedure  suggested  by  Liu  et  al.[5],  with  some  modifications,  was  used  to  prepare  the 
specimen  surface.  Figure  3.5  shows  schematically  the  steps  involved  in  this  process.  An 
optically  flat  glass  is  coated  by  vacuum  deposition  with  a  thin  silver  (or  aluminum)  film 
having  thickness  of  only  several  angstroms.  The  coated  glass  is  then  combined  with  the 
sample  using  an  epoxy  adhesive  to  glue  the  coated  surface  of  the  optically  flat  glass  to  the 
sample.  The  epoxy  adhesive  was  PC-1  Bipax  of  epoxy  resin  and  diethylentrianine  hardener 
obtained  from  the  Photoelastic  Division,  Measurements  Group  Inc.,  Raleigh,  N.C.  After  the 
epoxy  has  cured,  the  glass  is  peeled  off.  Because  the  bond  between  the  aluminum  film  and 
the  glass  is  weaker,  the  aluminum  coating  is  transferred  onto  the  sample  surface.  The  total 
thickness  of  the  epoxy  layer  and  the  coating  is  just  a  couple  of  microns.  Compared  with  the 
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sample  thickness,  this  layer  is  very  thin  and  will  not  affect  the  deformation  state  inside  the 
specimen. 


Figure  3.5:  Steps  involved  in  the  composite  surface  preparation. 

In  the  specimen  preparation,  two  points  are  remarkably  important.  First,  the  epoxy 
adhesive  should  be  free  of  viscoelastic  effects,  this  guarantees  the  flatness  of  the  surface 
minimizing  the  contraction  during  the  curing  time.  Second,  the  thickness  of  the  epoxy 
adhesive  should  be  minimized  to  avoid  shear  lag  effects  [13],  which  alters  the  measurement 
of  the  out-of-plane  deformations  in  regions  of  large  strain  gradients.  It  was  found  that 
slightly  polishing  the  specimen  surface  before  the  application  of  the  epoxy  adhesive  assists 
in  producing  a  uniform  epoxy  layer  giving  better  results. 

3.5  Application  of  the  CGS  Method  to  the  Quasi-Static 
Punch  Test 

First  the  quasi-static  punch  test  was  examined.  A  servohydraulic  testing  machine  was  used  as 
a  loading  frame.  The  specimen/punch  configuration  is  shown  in  figure  3.1.  The  materials  of 
the  specimens  were  PMMA  and  graphite-epoxy  composite.  Each  case  is  analyzed  separately 


below. 


3.5.1  Quasi-static  Punch  Problem.  Isotropic  Material 

Consider  a  rigid  punch  of  length  2a  acting  on  the  half  space  as  shown  in  figure  3.6(a). 
Neglecting  friction,  the  boundary  conditions  for  this  problem  are 

(ry{x,0)  =  0  for  |a;|  >  o, 

Txy{x,0)  =  0  for  |a;|<oo,  (3.24) 

v{x,  0)  =  Vo  for  |x|  <  a, 

in  addition  to  vanishing  stresses  at  infinity.  The  constant  displacement  of  the  punch,  no,  can 
be  arbitrarily  set  to  zero  with  the  proper  selection  of  the  origin  of  the  coordinates  x  —  y. 
The  total  normal  force  applied  by  the  punch  is  P. 


Figure  3.6:  (a)  Schematic  of  the  punch  problem  without  friction,  (b)  Two  external  semi¬ 
infinite  cracks. 

Note  that  the  boundary  conditions  (3.24)  for  the  rigid  punch  without  friction  are  the 
same  as  for  the  problem  of  two  external  semi-infinite  cracks,  figure  3.6(b),  loaded  symmetri¬ 
cally  by  two  remote  loads  P,  hence  both  problems  are  described  by  the  same  mathematical 
solution.  Using  complex  potentials  it  can  be  shown  [14]  that  the  stress  distribution  for  the 
frictionless  punch  problem  is  given  by 

a,  +  ay  =  2{4>{z)  +  W)}  =  4Re{</>(z)},  (3-25) 
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with 


(3.26) 


A.(  ^  ^ 

27r  Vi2-r^’ 

and  z  =  X  +  iy.  By  an  asymptotic  analysis  it  can  be  shown  that  the  stress  field  given  by 
equations  (3.25)  and  (3.26)  is  square  root  singular  near  the  punch  corners  (or  crack  tips) 
X  =  ±o,  and  a  stress  intensity  factor  can  be  derived  as 

.  K,  =  ^.  (3.27) 

VTTO 


Taking  the  derivative  with  respect  to  x  of  equation  (3.25)  gives 

=  -iKnp 
dx  V  TT 


(3.28) 


(^2  _  a2)3/2  (^2  _  a2)3/2 

where  Kj  =  P/y/im  has  been  used. 

Now  the  equation  governing  the  CGS  fringes  for  the  quasi-static  punch  test  may  be 
derived.  For  the  transmission  mode,  combining  equations  (3.20)  and  (3.28)  gives 


-iKj 


z 

-  a2)3/2 


(f2-a2)3/2j  cAA’ 


pm 


m  =  0,  ±1,  ±2, . . . 


(3.29) 


Experimental  Observations 

An  isotropic  PMMA  specimen  was  loaded  in  the  testing  machine  as  shown  in  figure  3.1.  A 
plate  of  PMMA  with  dimensions  203.2x50.8x6.35mm  was  used.  The  half-length  of  the  punch 
was  a  =  12.1mm.  A  50mm-diameter  collimated  He-Ne  laser  beam  (wave  length=632nm) 
was  employed  for  illumination.  Two  high-density  gratings  were  positioned  in  the  front  of 
the  specimen  with  the  gratings  parallel  to  the  a;— axis.  Therefore  information  regarding  the 
quantity  d{ax  +  o^jdx  was  obtained.  The  grating  pitch  was  p  —  0.0254mm  and  the  distance 
between  gratings  was  A  =  26mm.  From  previous  experiments  it  was  shown  that  the  optic 
constant  is  c  =lxl0“‘‘mm^/N  for  PMMA  [2]. 

Figure  3.7(a)  shows  the  expected  fringe  patterns  described  by  equation  (3.29)  for  a  given 
P.  Figure  3.7(b)  shows  a  photograph  of  a  real  interferogram  of  the  CGS  fringes  for  the  punch 
test.  Note  the  excellent  agreement  between  the  numerical  predictions  and  the  photograph. 
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(a)  (b) 


Figure  3.7:  (a)  Numerical  prediction  of  CGS  fringes  (constant  +  ay)/dx  values),  (b) 
Photograph  of  a  real  interferogram  of  the  CGS  fringes.  Isotropic  material,  transmission 
mode. 


CGS  Fringe  Interpretation 


To  determine  the  stress  intensity  factor  near  the  punch  corner,  the  photograph  of  the  CGS 
patterns  was  first  digitized  to  obtain  the  fringe  order  (m)  and  fringe  location  at  several 
points  (xfc,  yk)-  Then,  a  least  squares  scheme  was  used  to  fit  the  experimental  data  from  the 
photographs  to  expression  (3.29)  and  obtain  the  parameter  Kj. 

The  least  squares  method,  was  developed  as  follows.  First,  define  the  functions 


Fk  =  F{xk,yk)  = 


fl 

Zk 

/  TT 

{zk  — 

{Zk^  -  o2)3/2 

3/2 

fk 


and 


PWZfc  3/2 


where  Zk  =  Xk-^  m  and  r*  =  \zk\.  Then,  an  attempt  is  made  to  fit  the  expected  fringe 
pattern  (governed  by  equation  (3.29))  to  the  digitized  data  points  by  minimizing  the  error 
function 


m,)  =  E  [FkK,  -  , 


fc=l 


(3.30) 


where  N  is  the  total  number  of  points.  Note  that  since  there  is  only  one  parameter  to  be 
determined,  i.e.  Kj,  this  process  is  a  very  simple  least  squares  fit.  Thus,  employing  the 
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'  usual  methods  [15]  it  is  found  th^ 

:S  =  (3-31) 

It  is  noted  that  the  function  rp  heights  the  outer  lobes  of  the  interferograms  more  strongly 
than  the  inner  lobes. 

To  assure  the  accuracy  of  using  the  optical  CGS  technique  to  obtain  the  stress  intensity 
factor  near  the  punch  corner,  the  experimental  CGS  measurement  of  Kj  was  compared  with 
;  the  independently  measured  using  the  measured  load  and  equation  (3.27).  The  result  is 
illustrated  in  figure  3.8(a).  The  agreement  between  the  analytical  solution  and  the  optical 
experimental  measurement  is  excellent.  This  agreement  is  expected  since  the  theoretical 
lobes  predicted  by  equation  (3.29)  are  very  close  to  the  experimental  CGS  fringes  as  noted 
in  figure  3.8(b). 


Figure  3.8:  (a)  Comparison  of  the  experimental  measurements  of  Kj  with  the  analytical 
solution  for  the  punch  test  with  isotropic  materials,  (b)  Comparison  between  predicted  and 
observed  CGS  fringes. 


Uncertainty  Analysis 

Several  sources  of  uncertainty  are  involved  in  the  experimental  determination  of  the  stress 
intensity  factor  Kr.  However  as  will  be  shown  below  some  of  them  are  very  small  compared 
to  others  and  can  be  neglected. 
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Prom  standard  uncertainty  analysis  [15],  it  is  known  that  if  q{x, ...  ,2)  is  a  function  of 
several  variables,  then  the  uncertainty  of  q  in  terms  of  the  uncertainty  of  the  variables  is 


given  by  the  error  propagation  law 


(3.32) 


for  independent  random  errors. 


The  estimate  for  Kj  given  by  equation  (3.31)  is  a  function  of  the  measured  numbers 
ajv.  Thus  the  uncertainty  in  Kj  is  given  by  simple  uncertainty  propagation  in  terms 


of  the  uncertainties  in  0:1, ... ,  that  is 


-5a\  . 


(3.33) 


The  first  component  is  associated  with  the  fitting  process  error  and  the  second  with  the 
uncertainty  in  measuring  a.  It  can  be  shown  that  the  first  component  may  be  expressed  as 


^(dKjV  , 


EN  jp2  ’ 

k=l 


where 


N-1 


(3.34) 


(3.35) 


is  the  uncertainty  in  the  measurements  cki,  . . . ,  In  addition,  note  that 

dK,  SLi  SiiLi  -  2  ELi  FkT- 


where 


—  [(z2-a2)3/2  (42  _  ^2)5/2 J  * 


(3.36) 


(3.37) 


An  estimate  of  the  uncertainty  6K1  was  computed  according  to  equation  (3.33),  and  it 
was  found  that  the  second  component,  corresponding  to  the  uncertainty  of  a,  was  negligible 
in  comparison  with  the  first  one,  corresponding  to  the  fitting  analysis.  This  is  due  to  the 
fact  Sa  is  much  smaller  than  the  discrepancy  between  the  fitted  and  observed  fringe  pattern. 
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3.5.2  Quasi-Static  Punch  Problem.  Orthotropic  Material 

Consider  a  rigid  punch  of  width  2a  acting  on  a  half-plane  of  orthotropic  materials.  The 
mathematical  statement  of  the  problem  is  similar  to  that  of  the  last  subsection  for  the 
isotropic  case  (the  boundary  conditions  are  the  same). 

It  can  be  shown  [16]  that  the  in-plane  stress  components  for  this  case  are 

<7x  =  2Re{sl^{zi)  +  sl<f>{z2)}  (3.38) 

Cy  =  2Re  {^(^i) -f- ^(22)}  (3.39) 


where 

<l>{zi)  =  <t>{z2)  =  —^5(2:2), 

S2  —  Si  S2  —  Si 

and 

with  zi  =  x  +  siv,  =  S2y  and  5^,  (A;  ^  1, 2),  being  the  roots  of 


dllS^  +  (2ai2  +  CLqq)s^  +  022  —  O' 

The  coefficients  aij  are  related  to  the  compliance  matrix  Sij  by 

dij  —  Sij 


(3.40) 


(3.41) 


(3.42) 


for  plane-stress,  and  by 

=  (i,i  =  l,2.4,5.6) 

^33 

for  plane-strain.  The  compliance  matrix  is  given,  in  terms  of  engineering  constants,  by  [17] 


.Sn  = 

1 

Ei^ 

-  El  ’ 

544  = 

1 

IJ‘2z' 

S22  = 

1 

*^13  El  ’ 

^55  = 

1 

IJ-n 

S33  = 

1 

^3’ 

r,  _  “*^23 

<^23  —  rp  5 
£/2 

566  = 

1 

/^12 
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Prom  the  Generalized  Hooke’s  law  for  orthotropic  materials  it  can  be  shown  that 

(■z  —  S\z<^x  S2z(Ty  +  Szz(Tz'  (3.43) 


For  a  thin  plate  specimen,  it  is  assumed  that  plane-stress  deformation  conditions  prevail  and 
Gz  =  0.  Also,  according  to  the  generalized  plane-stress  assumption,  Cz  is  assumed  uniform 
through  the  specimen  thickness  and  the  out-of-plane  surface  displacement  w{x,y')  can  be 
related  to  the  in-plane  stress  components  through 

1  /•^/2  h 

w{x,y)  =  n  =  i:  +  SzzOy) .  (3-44) 

Z  J  ’-hjT.  ^ 

Taking  the  derivative  of  this  equation  with  respect  to  x  and  in  light  of  (3.38)  and  (3.39)  we 


h  fa. 


S2  Zi{Si3Si  +  S2z)  Si  Z2(SizS2  +  <^23) 

(a2-^2)3/2 


(3.45) 


where  Kj  =  P/y/Tra  has  been  used.  Finally  the  equation  governing  the  CGS  fringes  for  the 
quasi-static  punch  test  in  orthotropic  materials  may  be  derived.  For  the  reflection  mode, 
combining  (3.22)  and  (3.45)  gives 


S2  -3i(*S'i3Si  -h  <S'23)  Si  ^2(‘5'l3S2  +  Szz)  \  _  '^P 
S1-S2  (a2-z?)3/2  "si-S2  (a2-zi)3/2  J  “  2A 


(3.46) 


where  m  =  0,  ±1,  ±2, . . ..  Figure  3.9  shows  the  expected  CGS  fringe  patterns  described  by 
this  equation  for  a  given  Kj  and  the  material  properties  given  in  table  3.1. 


Experimental  observations  and  fringe  interpretation 


The  procedure  to  experimentally  examine  the  composite  material  was  similar  to  that  used 
in  the  isotropic  punch  test.  The  dimensions  of  the  specimen  were  100x50x6.3  mm.  The 
half-width  of  the  punch  was  a  =  15  mm.  The  x  gradient  was  considered,  i.e.,  dw/dx  was 
measured.  The  distance  between  gratings  was  A  =  44mm.  Figure  3.10  shows  a  photograph 
of  the  fringe  pattern  observed  in  the  composite  with  fibers  parallel  to  the  a;— axis.  Note  the 
similarity  between  the  numerical  predictions,  figure  3.9,  and  the  fringes  observed,  figure  3.10. 
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PunchProblem,  Orthotropic  material,  CGS  In  reflection 


X  (mm) 


Figure  3.9:  Numerical  prediction  of  the  CGS  fringes,  dw/dx,  for  the  punch  test  in  composite 
materials  with  fibers  parallel  to  the  x— axis. 


Figure  3.10:  Photograph  of  a  real  interferogram  of  the  CGS  fringes.  Composite  material 
with  fibers  parallel  to  the  x— axis. 

As  in  the  isotropic  case,  the  least  square  scheme  was  used  to  fit  the  experimental  data 
from  the  photographs  to  expression  (3.46)  to  find  the  parameter  Kj.  The  results  are  shown 
in  figure  3.11,  where  a  good  comparison  between  the  experimental  measurements  and  the 
analytical  expression  for  Ki  is  observed. 

An  estimate  of  the  uncertainty  SKj  was  computed  for  the  orthotropic  case.  From  the 
error  propagation  law  and  neglecting  the  uncertainty  5a  we  have 


5Ki  = 


(3.47) 


Quasl-statte  Punch  Test.  Composite  material 


Figure  3.11:  Comparison  of  the  CGS  experimental  measurements  with  the  analytical  solution 
for  the  punch  test  in  composite  materials. 

The  first  component  is  given  by  (3.34)  but  with  the  proper  definition  of  Fk  and  cxk-  The 
second  and  third  components  are  much  more  complicated  though.  Note,  for  example,  that 

dJQ  ^  ELi  Fl  ff -  2  ELi  Fkak  Fk^  ^3  ^3^ 

depends  on  dFk/dsi  and  this  depends  on  5si  which  depends  on  the  uncertainty  of  the 
material  properties  Ei,  E2,  //12  and  1^12.  Similarly  dFk/ds^  and  6s2  depend  on  Ei,  E2, 
and  U12. 

The  three  components  of  the  uncertainty  in  (3.47)  were  calculated.  However  it  was  found 
that  the  second  and  third  components  are  two  orders  of  magnitude  smaller  than  the  first  one. 
The  first  component,  due  to  the  fitting  process,  is  the  dominant  factor  in  the  uncertainty. 

3.6  Dynamic  Punch  Test 

The  dynamic  punch  test  was  conducted  on  PMMA  and  graphite  epoxy  composites.  The 
specimen  was  impacted  on  one  side  by  a  rigid  projectile  made  of  hardened  steel  as  shown 
schematically  on  figure  3.12. 
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3.6.1  Isotropic  Materials 

Experimental  Procedure 

The  stress  intensity  factor  for  dynamic  tests  was  examined  experimentally.  The  specimen 
geometry  is  shown  in  figure  3.12.  Specimens  are  made  of  PMMA;  the  mechanical  properties 
are  given  in  table  3.3.  Impact  of  the  specimen  is  achieved  by  a  cylindrical  steel  projectile 
shot  by  an  air  gun  as  shown  schematically  in  figure  3.13.  The  apparatus  consists  primarily 
of  an  air  gun,  a  Cordin  330  high  speed  camera  and  a  coherent  argon-ion  laser  pulsed  for 
10ns.  Once  the  projectile  is  launched  it  travels  along  the  barrel  and  activates  two  infrared 
detectors  placed  on  the  barrel  close  to  the  specimen.  The  signal  from  the  detectors  triggers 
the  camera  and  laser  controllers.  A  combination  of  delay  generators  make  it  possible  to  take 
pictures  at  the  exact  impact  time.  A  pulse  counter  limits  the  laser  to  only  80  pulses  of  10  ns 
each,  giving  80  images  per  test.  The  velocity  of  the  projectile,  Vp,  is  determined  by  recording 
the  time  it  takes  to  cross  the  infrared  detectors.  The  speed  of  the  camera  was  set  to  5.2/is 
between  pictures  or  192,300  frames/sec. 


Figure  3.12:  Specimen  and  projectile  geometry  in  the  dynamic  punch  test. 

Three  tests  were  performed  using  different  specimen/punch  dimensions.  In  all  the  cases 
the  thickness  of  the  PMMA  plates  was  6.3mm,  the  grating  pitch  was  p  =  0.0254  mm,  the 
distance  between  gratings  was  A  =30  mm,  and  a  50  mm-diameter  collimated  laser  beam 
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Table  3.3:  Mechanical  properties  for  PMMA 


E  (GPa) 

3.39  ±0.05 

u 

0.3  ±0.01 

c  (mm'^/N) 

le-4 

c?-'  (m/s) 

1904  ±15 

p  (Kg/m-*) 

1190  ±10 

Figure  3.13:  Schematic  illustration  of  CGS  set-up  in  transmission  mode  for  use  in  the  dy¬ 
namic  punch  test. 

was  used.  The  first  experiment  was  conducted  on  a  plate  with  dimensions  (see  figure  3.12), 
li  =127  mm,  I2  =50.8  mm  and  a  projectile  with  length  L  =145  mm  and  2a  =25.4  mm.  In 
this  case  the  projectile  impacted  the  specimen  at  the  velocity  Vp  =16  m/s.  For  the  second 
test  h  =100  mm,  k  =30  mm,  L  =145  mm,  2a  =30  mm  and  Vp  =18.5  m/s.  For  the  third 
test  h  =177  mm,  k  =50  mm,  L  =304  mm,  2a  =30  mm  and  Vp  =16  m/s. 

Experimental  Results 

Figure  3.14  shows  a  series  of  some  CGS  interferograms  for  the  first  test.  Note  that  the  size 
of  the  lobes  increases  and  then  decreases  with  time.  Such  a  behavior  is  expected  for  the 
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stress  intensity  factor  as  well.  An  analysis  of  each  photograph  was  performed  by  means  of 
procedure  similar  to  that  used  in  the  quasi-static  case.  This  lead  to  the  determination  of  the 
time  evolution  of  the  stress  intensity  factor  Ki{t)  for  the  punch  test. 


Figure  3.14:  Sequence  of  CGS  interferograms  for  the  dynamic  punch  test.  IVansmission 
mode,  isotropic  materials. 

Figure  3.15  shows  a  good  agreement  between  experimental  CGS  fringe  pattern  and  the 
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numerical  reconstruction  from  the  one  term  fitting  procedure.  The  time  evolution  Ki{t)  is 
shown  in  figure  3.16  for  the  first  experiment.  The  normalization  factor  is  Ko  =  aoCiyj2a/cd, 
where  Cj  is  a  material  dependent  parameter  given  in  equation  (3.9)  and  oq  =  pc^Vp.  The 
mass  density  p  and  the  dilatational  wave  speed  Cd  for  PMMA  are  given  in  table  3.3.  The 
particular  feature  for  this  experiment  is  that  the  pulse  duration  (as  determined  by  the  length 
of  the  projectile)  is  less  than  the  time  needed  for  the  dilatational  waves  to  reflect  from  the 
specimen  boundaries  and  reach  the  punch  corners.  As  discussed  in  section  2,  the  dynamic 
stress  intensity  factor  for  the  semi-infinite  punch  under  impact  load  aoHit)  applied  during 
a  time  interval  0  <t  <t*  (equation  (3.10))  is 

Ki(t)  =  (JoCi  [ViH{t)  -  -  t*)]  (3.49) 

where  Ci  is  a  material  dependent  parameter.  For  the  isotropic  case  in  mode  I  loading  [11] 

Cl  =  4^V",(1  -  c2/c2)/7r.  (3.50) 

Note  in  equation  (3.49)  that  the  stress  intensity  factor  AT/ (t)  due  to  a  pulse  a^Hit)  of 
duration  t*  is  equivalent  to  considering  two  semi-infinite  pulses  with  the  second  applied  after 
the  time  t*  and  of  opposite  sign  making  Kt{i)  decrease.  In  the  experimental  situation,  after 
the  punch  impacts  the  specimen,  a  compressive  dilatational  wave  propagates  through  the 
specimen,  and  since  the  boundaries  are  stress  free,  the  reflected  wave  from  the  opposite  side 
is  tensile  such  that  when  it  arrives  at  the  punch  corriers  it  makes  A"/(t)  decrease.  The  decay 
of  Ki{t)  in  both  cases  is  compared  in  figures  3.16-3.18. 

A  plot  of  equation  (3.49)  is  shown  in  figure  3.16  by  the  solid  line.  The  pulse  duration 
is  t*  =  49.4/is  which  is  the  time  for  the  dilatational  wave  to  travel  twice  the  punch  length, 
the  normalized  time  is  Cdt* /2a  —  3.35.  Note  that  C(i  denotes  the  dilatational  wave  speed  of 
the  specimen  which  is  different  from  that  of  the  projectile  since  they  are  made  of  different 
materials.  A  good  agreement  is  observed  between  the  solid  line  and  the  experimental  results. 
The  increase  and  beginning  of  decrease  of  Kj{t)  is  noted  in  the  experimental  results.  The 
solid  line  assumes  a  semi-infinite  plate  solution,  i.e.,  it  does  not  consider  wave  reflections 
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Figure  3.15:  Comparison  between  experimental  CGS  interferogram  and  numerical  simulation 
reconstructed  from  one  term  analysis. 

from  boundaries.  It  should  be  noted  that  the  comparison  of  the  experimental  results  is  made 
against  the  semi-infinite  punch  solution  (or  semi-infinite  crack).  Although  in  the  experiments 
the  punch  width  is  finite,  it  has  been  observed  [6]  that  the  effect  of  the  scattered  waves  from 
the  punch  corners  is  minimum. 

Figure  3.17  shows  the  dynamic  stress  intensity  factor  for  the  second  experiment.  Note 
again  a  good  agreement  in  the  loading  zone  between  experimental  results  and  the  theoretical 
prediction  for  the  semi-infinite  punch.  After  the  normalized  time  Cdt/2a  =  2,  decay  occurs 
on  Kj{t)  due  to  the  arrival  at  the  punch  corners  of  the  reflected  waves  from  the  side  of  the 
specimen  opposite  to  the  side  impacted  by  the  punch. 

The  third  experiment  is  similar  to  the  second  test  in  the  sense  that  the  decay  of  Kj{t), 
at  normalized  time  Cdt/2a  =  3.33,  is  due  to  the  arrival  at  the  punch  corners  of  the  reflected 
waves  from  the  opposite  side  of  the  specimen,  figure  3.18.  Due  to  the  specimen/punch 
dimensions  for  this  experiment,  the  length  of  the  pulse  Cdt* /2a  =  5.93  is  long  enough  such 
that  it  does  not  affect  the  behavior  of  A/(t),  i.e.,  it  does  not  contribute  to  the  decay  of  Ki{t) 
as  in  the  first  experiment. 

The  uncertainty  analysis  for  this  case  requires  one  to  estimate  6  ,  that  is,  the  un- 
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Dynamic  SIF,  test  1 


Figure  3.16:  Normalized  dynamic  stress  intensity  factor  for  the  punch  test.  First  test,  short 
pulse  duration. 

certainty  of  the  normalization  factor  should  be  considered  in  the  expression 


5Ki  is  determined  as  in  the  quasi-static  punch  problem.  Neglecting  the  uncertainty  5a  we 
found  from  the  definition  of  Kq  that 


A  successive  application  of  the  error  propagation  law  is  needed  to  express  5cq,  SCj  and 
5cd  in  terms  of  the  uncertainty  of  the  material  properties.  After  the  proper  calculation  it 
was  found  that  that  6  (^)  is  different  from  ^  by  less  than  5%,  which  means  that  that 
the  uncertainty  in  the  normalization  factor  may  be  neglected.  Again,  the  largest,  dominant 
uncertainty  is  that  due  to  the  fit. 

3.6.2  Orthotropic  Materials 

Experimental  Procedure 

The  dynamic  punch  test  was  conducted  on  unidirectional  fiber-reinforced  composite  materials 
to  measure  the  time  evolution  of  the  stress  intensity  factor  Ki{t).  The  graphite-epoxy 
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Dynamic  SIF.  test  2 


Figure  3.17:  Normalized  dynamic  stress  intensity  factor  for  the  punch  test.  Second  test. 

composite  used  in  the  dynamic  experiments  is  the  same  as  that  employed  in  the  quasi-static 
test,  and  the  relevant  properties  are  given  in  tables  3.1  and  3.2. 

The  CGS  set-up  in  reflection  for  the  dynamic  punch  test  is  shown  schematically  in  flgure 
3.19.  Note  the  addition  of  the  beam  splitter  and  the  different  laser  beam  path  with  respect 
to  the  set-up  for  transmission  mode,  figure  3.13. 

With  reference  to  the  figure  3.12,  the  dimensions  of  the  specimen/projectile  used  in  the 
test  are  li  =  101.6mm,  I2  =  50.8mm,  h  =  6.5mm,  L  =  304mm,  2a  =  30mm.  The  fiber 
orientation  of  the  specimen  was  along  the  a:— axis.  The  x— gradient  was  considered,  i.e., 
dw/dx  was  measured.  The  projectile  impacted  the  specimen  at  the  velocity  Vp  =  11.5m/s. 
The  grating  pitch  was  p  =  0.0254mm  and  the  distance  between  gratings  was  A  =  36mm. 
The  speed  of  the  camera  was  set  to  4/iS  between  pictures  or  250,000  frames/sec. 

Experimental  Results 

A  sequence  of  pictures  displaying  the  CGS  fringe  pattern  for  the  dynamic  punch  test  in 
composites  is  shown  in  figure  3.20.  An  analysis  of  each  photograph  was  performed  by  means 
of  a  similar  procedure  to  that  used  in  the  quasi-static  case  for  composites,  and  thus,  the 
stress  intensity  factor  time  evolution  Ki{t)  is  determined. 
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Dynamic  SIF.  test  3 


Figure  3.18:  Normalized  dynamic  stress  intensity  factor  for  the  punch  test.  Third  test. 


A  comparison  between  the  experimental  CGS  fringes  and  the  numerical  reconstruction 
from  the  one  term  fitting  procedure  is  shown  in  figure  3.21,  note  the  good  agreement  in  both 
sets  of  firinges.  The  normalized  dynamic  stress  intensity  factor  for  the  punch  test  is  shown  in 
figure  3.22.  The  normalization  factor  is  Kq  =  where  oq  =  pCd,yVp,  being  Cd,y 

the  dilatational  wave  speed  for  wave  propagation  along  the  ^-axis  and  is  given  in  table  3.2. 
The  normalized  time  is  Cd,yt/2a.  An  increase  in  Ki{t)  is  noted  and  then  a  decrease  occurs 
when  the  dilatational  wave  reflected  from  the  opposite  side  of  the  specimen  arrives  at  the 
punch  corners.  The  solid  line  corresponds  to  the  stress  intensity  factor  for  the  semi-infinite 
punch  problem  under  impact  load  applied  during  a  time  interval  0  <t  <t*  (equation 

(3.10)) 

Kj{t)  =  aoCi  [ViH{t)  -  y/t^H{t  -  f)]  (3.51) 


where  Ci  is  a  material  dependent  parameter  derived  by  Rubio-Gonzalez  and  Mason  [10]  and 
given  in  equation  (3.8), 


Cl  =  2, 


2c,e 


(3.52) 


The  uncertainty  <5  was  estimated  following  the  procedure  described  for  the  dynamic 
punch  test  in  isotropic  materials.  Again  it  was  found  that  is  not  very  different  from 
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Figure  3.19:  Schematic  illustration  of  CGS  set-up  in  reflection  mode  for  use  in  the  dynamic 
punch  test  of  composites. 

i.e.,  the  uncertainty  in  the  normalization  factor  can  be  neglected. 

3.7  Conclusions 

Like  in  dynamic  fracture  mechanics,  the  dynamic  stress  intensity  factor  Ki  has  been  intro¬ 
duced  to  characterize  the  severity  of  the  singular  stresses  around  the  punch  corners  in  the 
punch  test.  Quasi-static  and  dynamic  punch  tests  were  performed  on  PMMA  and  unidirec¬ 
tional  graphite-epoxy  composites  to  measure  Kj.  Laser  interferometry  was  used  to  measure 
in-plane  stresses  (transmission  mode)  and  out-of-plane  displacements  (reflection  mode) ,  and 
then  estimate  the  stress  intensity  factor.  In  the  dynamic  case,  a  high  speed  photography 
technique  was  employed  to  capture  the  transient  response  of  the  specimen  and  measure  K{t) 
just  after  the  impact.  Each  photograph  of  the  CGS  fringe  pattern  was  digitized.  A  fitting 
procedure  based  on  the  least  squares  method  was  used  to  estimate  Kj  from  the  fringe  pat¬ 
tern.  In  all  the  cases  a  good  agreement  between  the  measurements  of  Kj  and  the  theoretical 
predictions  was  found.  The  theoretical  prediction  for  the  dynamic  case  on  composites  is  just 
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t=36  ^iS  I 


Figure  3.20:  Sequence  of  CGS  interferograms  in  reflection  for  the  dynamic  punch  test  in 
composites. 

an  extension  of  the  solutions  developed  for  orthotropic  materials.  [6,  10] 

An  estimate  of  the  uncertainty  in  the  experimental  determination  of  Kj  was  made.  It 
was  found  that  the  dominant  contribution  to  the  uncertainty  5Ki  was  that  associated  with 
the  fitting  procedure. 
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Figure  3.21:  Experimental  CGS  interferograms  and  numerical  simulation  reconstructed  from 
one  term  fitting  analysis.  Composite  materials,  (a)  At  t  —  12fis  and  (b)  at  t  =  20/iS. 


Dynamic  SIF,  punch  test,  composite  material 


Figure  3.22:  Normalized  dynamic  stress  intensity  factor  for  the  punch  test.  Composite 
materials 
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